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ABSTRACT 


The electrical behaviour of microwave varactors under both 
steady-state and large-signal transient conditions is investigated. 
Numerical iterative methods are used to solve the one-dimensional two- 
carrier transport aan which are taken to describe the behaviour of 
these devices under both steady-state and transient conditions. The 
solutions for the steady~state and time-dependent cases are obtained by 
modifying and extending numerical methods, already available in the 
literature, to cover the present cases. Considerable attention is given 
to the numerical analysis of the steady-state and time-dependent cases in 
order to achieve a numerical algorithm sufficiently sound to cope with 
several severe difficulties of the problem. Two diodes, taken to 
represent two distinct kinds of microwave varactor, are considered. 
Results have been computed for four external excitations, namely: steady- 
state voltage excitation, "turn-on", "turn-off" and sinusoidal excitation 
in a broadband circuit. A slight charge accumulation in the vicinity of 
ian? wi space charge region during the "turn-off" transient in the 
>’ =S Oe fe nt step-recovery diode, and an inductive behaviour of these 
devices under large-signal sinusoidal excitation, are two of the many 
interesting aspects revealed by the results. The microwave measurement 


techniques used for the varactor measurements made are also presented and 


discussed. 
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Spatial distribution of electrostatic potential 
and electron and hole quasi~Fermi levels for 
diode No. 2 under steady-state conditions 
FOLreVy = 0.6 VOLT 


Spatial distribution of electrostatic potential 
and electron and hole quasi-Fermi levels for 
diode No. 2 under steady-state conditions 

for Va =) 10 VOLE 


Spatial distribution of net recombination rate 
for diode No. 2 under steady-state conditions 
for v= 025 -voLt 


Spatial distribution of net recombination rate 
for diode No. 2 under steady-state conditions 
for Vay = 0.8 volt 


Spatial distribution of net recombination rate 
for diode No. 2 under steady-state conditions 
for Va = 1.0 volt 
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Spatial distribution of carrier mobilities for 
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Spatial distribution of carrier mobilities for 
diode No. 2 under steady-state conditions 
for Va = 0.5 volt 


Spatial distribution of carrier mobilities for 
diode No. 2 under steady-state conditions 
for Va) =.0.8 volt 


Spatial distribution of carrier mobilities for 
diode No. 2 under steady-state conditions 
LOPON eC 1sO8 volt 


Spatial distribution of electron and hole 
currents for diode No. 2 under steady-state 
CondmedonssforeV, = 0.5 volt 


Spatial distribution of electron and hole 
currents for diode No. 2 under steady~state 
conditions for Vy, = 0.8 volt 


Spatial distribution of electron and hole 
currents for diode No. 2 under steady-state 
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Terminalncurrent of diode No. 2 as’ a ftnction 
of terminal voltage (steady-state conditions) 


Incremental capacitance of diode No. 2 as a 
function of terminal voltage (steady-state 
conditions) 


Turn-on voltage response of diode No. 1 for 
a laree, current step 


Turn-on voltage response of diode No. 1 for 
the initial phase of the turn-on transient 


Evolution of the electric field distribution 
in diode No. 1 during the turn-on transient 
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Evolution of the electron density distribution 
in diode No. 1 during the turn-on transient 
(values of time instants indicated in 
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Eyolution of the hole density distribution 
in diode No. 1 during the turn-on transient 
(values of time instants indicated in 
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Turn-on voltage response of diode No. 2 for 
a large current step 


Turn-on voltage response of diode No. 2 for 
the initial part of the turn-on transient 


Evolution of the electron density distribution 
in diode No. 2 during the turn-on transient 


Evolution of the hole density distribution in 
diode No, 2 during the turn-on transient 


Evolution of the electric field distribution 
in -diode No. 2 during the initial part of the 
turn-on transient 


Evolution tof the electric field distribution 
in diode No, 2 during the later part, of, the 
turn-on transient 


Evolution of the displacement current 
distribution in diode No. 2 during the initial 
part of the turn-on transient 


Evolution of the displacement current 
distribution in diode No. 2 during the later 
part of the turn-on transient 


Evolution of the electron drift current 
distribution in n'-region of diode No. 2 
during the turn-on transient 


Circuit model used for the study of the large~- 
signal turn-off response of the diodes 


Switching waveform applied to the diodes 


Terminal voltage as a function of time during 
the initial period of the reverse transient, 
for diode No. l 


Evolution of the electric field in the vicinity 


of the junction during the initial period of 
the turn-off transient, for diode No. l 


Evolution of the electric field throughout 
diode No. 1 during the initial period of the 
reverse transient 


Spatial distribution of hole and electron 
currents for a particular time instant 
(= 0255. X 1073) during the reverse transient, 


for diode No. l 


Diode terminal voltage and current as a 
function of time during the reverse transient, 
for diode No. l 


Hole density distribution at various instants 


of time during the reverse transient, for 
diode No. 1 
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Electron density distribution at various 
instants of time during the reverse transient, 
Lor -diode, No. 1 


Electron density distribution at various 
instants of time during the reverse transient, 
for diode No. 1 


Terminal voltage and current as a function of 
time during the initial period of the reverse 
transient, for diode No. 2 


Electric field distribution at various instants 
of time (values indicated in Fig. 9.9) during 
the initial period of the reverse transient, 
for diode No. 2 


Terminal voltage and current as a function of 
time during the reverse transient, for 
diode No. 2 


Electric field distribution at various instants 
of time (values indicated in Fig. 9.11) during 
the reverse transient, for diode No. 2 


Hole density distribution at various instants 
of time (values indicated in Fig. 9.11) during 
the reverse transient, for diode No. 2 


Electron density distribution at various 
instants of time (values indicated in Fig. 9.11) 
during the reverse transient, for diode No. 2 


Spatial distribution of hole and electron 
currents for a particular instant of time 

(= 0.29 X 10-4) during the reverse transient, 
for diode No, 2 . 


Circuit model used for the study of the large- 
signal, sinusoidal excitation response of 
the diodes 


Applied waveform in the circuit of Fig. (10.1) 


Terminal voltage, applied voltage and terminal 
current versus time for the first five cycles 
after the signal is applied to diode No. 1 


Spatial distribution of hole and electron 
densities and of electric field at time instant 
ty (value indicated in Fig. 10.3) for diode 
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Spatial distribution of hole and electron 
densities and of electric field at time instant 
to (value indicated in Fig. 10.3) for diode 
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Spatial distribution of hole and electron 
densities and of electric field at time instant 
te (valueyindicated in Fie. 10.3) for dicde 
NOww. 


Spatial distribution of hole and electron 
densities and of electric field at time instant 
tin(value indicated in Fie. 10.3)) for, dicde 

No. l 


Spatial distribution of hole and electron 
densities at various instants of time during 


the fourth cycle (values indicated in Fig. 10.3), 


FOr artode NO. ol 


Terminal voltage and current versus time for 
about the first cycle of the applied signal 
for diode No. 2 


Spatial distribution of hole and electron 
densities.and of electric’ field at time instant 
tj (value indicated in Fig. 10.6) for diode No. 


Spatial distribution of hole and electron 
densities and of electric field at time instant 
t2 (value indicated in Fig. 10.6) for diode No. 


Spatial distribution of hole and electron 
densities and of electric field at time instant 
t3 (value indicated in Fig. 10.6) for diode No. 


Spatial distribution of hole and electron 
densities and of electric field at time instant 
ty, Walue Indicated in. Fie. 10.6) for diode No. 


Spatial distribution of Nolte endvelectron 
densities and of electric field at time instant 
t5 (value indicated in Fig. 10.6) for diode No. 


Spatial distribution ‘of hole and electron 
densities and of electric field at time instant 
t6 (value indicated in Fig. 10.6) for diode No. 


Terminal voltage and current versus time for 
the initial transient period for diode No, 2 


Sectional view of junction between a rectangular 


waveguide (TEjg mode) and coaxial line 
Equivalent case of configuration shown in 

Pio Cla, 1) 

Equivalent circuit for the configuration shown 
inj Fig. Cll2Z) 


Complete equivalent circuit for waveguide 
shunted by a small-diameter dielectric post and 


an impedance Z, presented by a coaxial line 
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Simplified equivalent circuit applicable for 
a centrally-mounted diode in the waveguide 


Cross-section of a typical packaged 
varactor diode 


The reduced-height rectangular waveguide 
diode mount 


Circuit arrangement used for measurements on 
waveguide mounted varactor diodes 


Slotted-line calibration curve at 6.757 GHz 


Packaged varactor diode equivalent circuit 
at microwave frequencies 


Simplified equivalent circuit of packaged 
varactor diodes at lower microwave frequencies 


Low-frequency terminal capacitance versus 
bias for the step-recovery diode measured 
(HP. No. 5082-0251) 


Step-recovery diode (HP. No, 5082-0251) _ 
impedance for various bias voltages derived 
by Harrison's method at 6 GHz 


Step-recovery diode (HP. No. 5082-0251) 
impedance for various bias voltages derived 
by Harrison's method at 7 GHz 


Calculated incremental capacitance of diode 
No. 2 and low frequency terminal capacitance 
of step-recovery diode (HP. No. 5082-0251) 
as functions of forward-bias voltage 
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CHAPTER I 
INTRODUCTION 


The varactor diode is a semiconductor device that is used as a 
variable-reactance circuit element. The variable reactance is provided 
by the junction capacitance, which is a nonlinear function of applied 
Peseasal es The varactor is therefore a nonlinear element, and this 
nonlinearity can produce substantially different effects. One of these 
is the switching or modulation of a microwave signal through variation 
of the reactance by means of an externally applied voltage. In a second 
usage, the nonlinearity causes the generation of harmonics of an applied 
microwave signal. In a third case, two microwave signals of different 
frequency may be applied, resulting in parametric amplification, 
up-conversion or down-conversion of one of the signals. 

The purpose of the present work is to investigate both the 
Static and dynamic behaviour of microwave varactors. The popular use of 
overdriven Phe rar er ce necessitated the need of an analysis which covers 
forward-biased conditions in a relatively general form. Moreover, 
the response of a varactor driven by ideal steps of current and voltage 
is of considerable help in understanding its dynamic behaviour. In this 
study, four basically different forms of external excitations are chosen 
to be analysed. First is the behaviour of a varactor when a forward- 
biasing steady-state terminal voltage is specified. Second, the external 
excitation chosen is a large current step applied to the varactor 
initially in thermal equilibrium condition ('turn-on' transient). The 


third case deals with the situation that occurs when a reverse~-biasing 
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potential is suddenly applied to a varactor initially in the steady-state 
forward-biased condition ('turn-off' transient). The fourth aspect that 
is selected to be analysed is the case of a reverse-biased varactor 
excited by a sinusoidal microwave signal (frequency in the GHz range) 
that is large enough to drive the varactor into forward-biased condition 
for a part of the cycle. 

Basic concepts in the theory of p-n junctions were first 
presented by Bhoeialc ate Shockley's steady-state and small signal 
analysis was based upon the simplifying assumption that the space-charge 
layer of a p-n junction is completely depleted of mobile carriers. This 
assumption eliminates from Poisson's equation the distribution term for 
mobile holes and electrons and thereby reduces a difficult system of 
differential equations to a single, mathematically tractable differential 
equation. Despite this simplifying assumption, many conclusions derived 
from the depletion layer approximation are verified by laboratory 
experiments. 

Shockley's depletion layer approximation contains severe 
restrictive features rendering it inapplicable to several modes of 
junction operation. In particular, to assume that the space-charge layer 
is free of mobile charge carriers implies that such a theoretical 
approximation is not applicable to forward-biased p-n junctions, nor to 
reverse-biased junctions containing a large electric current. Although 
these limitations have been recognised by many workers, no rigorous 
analytical treatment of the problem has been reported. Instead, other 
ee eeness have been used to simplify the mathematical equations 
characterising a semiconductor junction; each approach, however, has 


introduced a new set of uncertainties concerning the analytical 
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limitations arising from the use of simplifying assumptions. 

A number of orca have studied the turn-on and turn-off 
transient in p-n junctions. These investigations have been conducted 
upon analytical models which neglected many physical mechanisms 
encountered in practical devices. Moreover, only highly asymmetrical p-n 
junctions were considered. nOnGE: Cena gave a second order analytical 
treatment for p-i-n rectifiers. This analysis also neglected many 
physical mechanisms which are encountered in actual practice. To our 
knowledge, no comprehensive analytical treatment has been reported for a 
p-n junction excited by a sinusoidal signal large enough to forward-bias 
it for part of the cycle. The presence of the passive elements 
surrounding the p-n junction devices in the actual circuitry further 
complicates the already formidable mathematical problem. 

Some of the assumptions generally adopted in almost all of the 
analytical work for the one-dimensional model (to which attention will 


be limited in this thesis) are the following: 


(a) Separation of the structure into regions with sharp 


boundaries, either fully depleted of mobile carriers, or 


space-charge neutral; 


(b) Postulation of explicit boundary conditions on relevant 
quantities in the interior of the device and at the 
interface between the depleted and neutral regions of 


assumption (a); limited choice of boundary conditions at 


the external contacts; 
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(c) Limitation of the doping profile to very special cases 
(mostly step and linear distributions) and to particular 
quantitative values (mostly highly asymmetric impurity 


distributions) ; 


(d) Simplification of the dependence of carrier mobilities upon 


the electric field, doping and scattering phenomena; 


(e) Limitation of carrier recombination laws to the low-level 


linear case. 


Numerical methods, with the aid of high-speed digital computers 
represent an alternate approach to the problem, the final aim being the 
achievement of a solution of the most general character without the 
conventional assumptions. Serious difficulties are also present in a 
numerical investigation, and have prevented most of the currently 
available numerical solutions from having the general character desired. 
The first comprehensive and general numerical procedure for the solution 
of the steady-state problem is presented by Garnet and is applied to 
the solution of transistors. nate later on exposed some of the 
limitations of Gummel's analytic formulation, slightly modified it and 

20 
applied it to very short junction structures. Kennedy et al. used 
numerical techniques to solve the zero-current model of the linearly- 
graded p-n junction. Later, ava reported a numerical method of 
solution of p-n junctions under arbitrary transient conditions and 
applied it to his simpler model having zero-recombination and constant 
mobilities to obtain solutions for the case of a highly asymmetrical step 


p-n junction. Another algorithm for time-dependent solutions of p-n 
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junction devices is reported by Scharfetter and Gomes they applied 
it to solve for the large-signal behaviour of IMPATT diodes. 

In this thesis numerical solutions for a relatively general 
physical model, considered sufficient for microwave varactors, are 
obtained to study the device behaviour under both static and dynamic 
conditions. The solutions for steady-state conditions are based on the 
algorithm furnished by cunmer The algorithm is slightly modified to 
make it more efficient at higher forward biases where Gummel's original 
algorithm tends to fail. The time-dependent solutions are based upon an 
aes of the algorithm furnished by DeMee for our more general 
physical model. Suitable numerical techniques are selected to deal with 
various critical situations inherent in the solution for our particular 
device. All the reported solutions are for one-dimensional cases. 
Microwave measurements of the static properties, particularly the series 
resistance, are also undertaken. 

The presentation is divided into eleven sections. Chapter II 
discusses various design considerations for microwave varactors and 
describes the particular device structures chosen to represent two 
distinct types. Chapter III describes the important physical mechanisms 
included in the physical model adopted and discusses implications of the 
various assumptions made. Chapter IV describes the fundamental equations 
and boundary conditions that determine mathematically both the steady- 
state and time-dependent problems. Chapter V describes the details of 
the analytical formulation and the numerical techniques adopted for the 
steady-state case. Chapter VI presents the analytical formulation and 
numerical techniques suitable for the achievement of time-dependent 


solutions for the present physical model. It is an extension of the 
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numerical algorithm presented by peter” for his simpler physical model. 
This chapter also discusses various sources of numerical error present 

in the algorithm and various trade-offs between different errors to 
decrease the overall error without increasing computation time. 

Chapter VII presents the results for the steady-state case: terminal 
properties and quantities in the interior of the device are displayed; 
the results are briefly discussed. Chapter VIII presents the results for 
the case when a large step of current is applied to the varactor initially 
in thermal equilibrium (turn-on transient). A brief discussion on the 
results is also offered, Chapter IX presents the results for the case 
when a Pee hie potential is applied to the varactor initially in 
the forward-biased condition (turn-off transient). The results are also 
discussed. Chapter X presents the internal and external behaviour of the 
Varactor in a broad-band multiplier circuit. Chapter XI describes the 
various microwave techniques employed to measure the static character- 
istics:.of.a varactor. Conclusions are given in Chapter XII; some of the 
several immediate applications for the analysis under more general 
conditions and the possibilities of extensions to more complex situations 


are briefly discussed. 
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CHAPTER II 
MICROWAVE VARACTOR DIODES 


This chapter describes the choice of physical structures 
representing the main categories of microwave varactor diodes. A simple 
mathematical model is formulated to predict the impurity distribution 
inside the diode structure. The parameters in the mathematical model 
are then chosen so that the desired impurity distribution is obtained. 
At the end of this chapter, the impurity distributions inside the diodes 
are displayed and various structural parameters of the diodes are 


bapuLtated. 


2, Lilatroduction 

Microwave varactor diodes can be classified into two main 
categories. The diodes belonging to the first category have a varying 
Bei tarine in the reverse voltage direction. The movement of the 
potential wall in the depletion region uncovers charge continuously, and 
results in the capacitance change. This capacitance change in the 
reverse direction is used in varactor multipliers to obtain harmonic 
currents, which are related to the order of non-linearity in this 
Capacitance, This class of microwave varactors will be mentioned as 
‘varactor diodes' in this thesis. The diodes belonging to the second 
category exhibit a relatively constant reverse-bias capacitance and a 
very large diffusion capacitance in the forward direction, resulting in a 
much greater variation of the total capacitance of a varactor. This 


~ * ' ue * ' 
class of microwave varactors will be mentioned as ‘step-recovery diodes 


io chis: thesis. 
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2.2 Design Considerations 

Various useful varactor diode structures are described by Chang 
and eee They show that hyper-abrupt junctions, and "Pagoda" 
structures exhibit a relatively large capacitance variation. However, 
when compared to a step~junction having the same breakdown voltage, the 
improvement due to increased capacitance-voltage sensitivity is offset, 
or more than offset by the concomitant increase in the R-C product. 
Therefore a heavily doped step-junction can be taken to represent micro- 
wave varactors. It is almost impossible with the presently developed 
technology to construct a p? - ae diode with perfectly abrupt steps in 
the impurity concentration at the boundaries, Hence a doubly-diffused 
epitaxial diode in which the diffusion is very deep, i.e., the junction 
is formed by the impurities diffused in from the surface and substrate, 
is chosen as a typical microwave varactor. The impurity concentration at 
the substrate and at the surface are assumed to be very high to ensure 
low series resistance. 

In the step-recovery diodes a constant reverse-bias capacitance 
and a very high diffusion capacitance are required. The diffusion 
capacitance is produced by the injection of minority carriers during the 
forward bias excursion of the varactor pump voltage and the subsequent 
withdrawal of this charge during the reverse bias portion of the waveform. 
In a step-recovery diode the following requirements must be met: Firstly, 
the junction must be capable of minority carrier injection, i.e.5ethe 
injection of holes from a p-type region into an n-type region, and/or 
electrons in the other direction. The second basic requirement is that 
the injected charge be recoverable during the reverse-bias half of the 


cycle. In order for this to occur, two conditions are necessary; first, 
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the carriers must not recombine appreciably before they are sainakeahane 
and second, the carriers must not diffuse so far from the junction that 
they cannot be retrieved. The first condition is easily met by most 
silicon epitaxial microwave diodes, since the typical minority-carrier 

iS : : . =o -6 
lifetime in such a material ranges from 10 to > X 10)» seconds. 

; : : nee 2) 

Moreover, the diode can be subjected to higher lifetime treatment to 
achieve still higher minority-carrier lifetime. Ga-As will fulfil this 
condition at higher microwave frequencies because minority-carrier life- 

; : =o A 
time is 10 seconds or less. Nearly all step-recovery diodes used at 
microwave frequencies below millimeter range are made of silicon. To 
Satisfy the second of the above conditions, i.e., to prevent the minority- 
carriers from diffusing beyond the point of no return, large built-in 
electric fields on both sides of the junction are necessary. Also, 
because of the finite time required for the carriers to cross the 
depletion-region and current-limiting effects of the field produced by 
the injected carriers themselves, the depletion region must not be too 

: Z = e Sea : ; A 
wide. The ideal p - i-n diode is one with perfectly abrupt steps in 
impurity concentration at the boundaries with the i-layer, having a width 
between 0.5 and 4 microns. The epitaxial growth and diffusion occur in 
finite times at the diffusion temperature so that i-layer boundaries are 
graded. The device has to be fabricated such that the impurity profile 
is very steep. Moll et Sime have stated that if the width of the 

: ; : ESN ee 

gradient region measured to an impurity concentration of 10° cm is 
less than 1/4 to 1/2 of the i-layer width, there is insignificant 
rounding-off and relatively little dependence of transient behaviour on 


the charge injection level. This criterion has been utilised in choosing 
o 
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the impurity doping profile in this study. Undoped growth of silicon 
generally results in an n-type material. Thus the high-resistivity 
center region of the step-recovery diode analysed is assumed to be 

, 148 , ‘ : 
slightly donor doped (107° cm ~). The steep impurity profiles and the 
narrow junctions required for fast and complete recovery of the injected 
charge necessarily result in a lower breakdown voltage. The high- 
resistivity center region decreases the cut-off frequency of the device, 
- : . , 1 : : 
since the effective series resistance of a varactor is larger when its 
charge-storage capacitance is used than when operation is confined 
Strictly to’ dielectric capacitance. 
Another consideration which is important is the choice of 
junction area. We know that the junction capacitance and charge stored 
, : ets 24 2 ‘ , 
during forward-bias conditions are proportional to the junction area. 
; ; : , ; 24 
In step-recovery diodes the reverse-bias capacitance is undesirable  ; 
therefore, from this point of view, small junction area is desirable. As 
a result, a compromise has to be made between the charge-stored which 
decreases as the area is decreased and the reverse-bias junction 
capacitance which can be tolerated. The higher mobility of free 
electrons compared to that of the holes results in a lower resistivity 


for the same doping level. For this reason, the substrate in epitaxial 


mesa varactors is usually chosen to be n-type. 


2.3 The Fabrication Process and its Mathematical Model 


Poet taniication Process 


Fabrication seapeee include epitaxial growth of the i-n 
ms ‘ s + 
boundary and subsequent growth or diffusion of the i- p boundary. 


These steps are followed by more or less conventional semiconductor 
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fabrication steps to shape it into a mesa-type structure. The diffusion 
time and temperature are chosen so as to yield the required impurity 
gradient. In the case of varactor diodes, deep diffusion is required, 
while in step-recovery diodes shallow diffusion is needed. In some cases 
the process is simplified and diffusion of donor impurities from the 
substrate and acceptor impurities from the surface are carried out 
Meeetdecieny? Boron is one of the most useful p-type impurities and 
Hs chosen’ as the p-type impurity for our calculations. The substrate is 
taken to be phosphorous-doped. Phosphorous and boron have the same 
diffusion coefficient for 1000 - 1200°c, the most usual diffusion 
temperatures. For a varactor, where deep~diffusion is required, the 
substrate can be doped with arsenic (which has a diffusion coefficient 
whose value is 1/3 tia tor boron) to yield a deeper junction from the 
p-type impurity side. This in turn makes the ohmic contact to the 
surface easier to produce. 


2.3.2 The Mathematical Model 


28 : a 
Fick's first and second laws are taken to govern the diffusion 


phenomenon. 

f= ee) 

ox 

on pele &. 2) 

ot ox 
where 

fe= nartrcle Lux through a plane parallel to the 

—2 -1 
surface (cm sec ), 
a 2 —l 
D = diffusion coefficient (cm™ sec ), 
. —3 
N = impurity concentration (cm ~). 
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The first equation describes the rate of permeation of the diffusing 
impurities through a unit cross-sectional area of the host medium. The 
second equation specifies the rate of accumulation of diffusing 
impurities at different points in the host medium as a function of time, 
The in-diffusion from the surface can be considered as if the 
epitaxial layer were of infinite thickness and its background impurity 
concentration negligible compared to the surface concentration. Then, 


for boundary conditions: 


N(O,t) Eau No3 N(x,0) = 0 wiege 2 oD) 


where 


N. = impurity concentration at the surface (ones) and 
x = distance from the surface (cm). 
The Snbavion © of Eq. (2.2) for a semi-infinite medium and boundary 


conditions of Eq. (2.3) yields the p-type impurity concentration- 


distribution: 


NGot) = 5, (1 = ont ~ ee) 
a S 


The n-type impurity distribution due to out-diffusion from the 
substrate can similarly be considered to occur as if the substrate and 
the epitaxial layer were of semi-infinite width; the background impurity 
concentration can be considered to be negligible compared to the 
substrate impurity concentration. Then, solving Eq. (2,2) forthe 
boundary conditions: 


N(x,0) = No srs = ae 
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and 


N(x,0) = 0 O< x < + 


> 


the donor impurity distribution can be expressed as follows (using the 


Petein fOr x as in Eq; 2.4) 


No a 
N4 Ost) = 2(1 ~ aa »o0 (2.5) 
where 
os epitaxial layer thickness, 
No = impurity concentration of the substrate Gat and 


Di» (D.) = diffusion constant of p-type (n-type) impurity. 


2 
Detailed derivation of Eqs. (2.4) and (2.5) is given by Boltak oe The 
net impurity distribution in the epitaxial layer after diffusion can be 
obtained from Eqs. (2.4) and (2.5) and its original impurity 


distribution, Thus: 


N x, - x 
NG) = -a,(1 - ext z= yeas +N, Sn as) 
" avd. t ‘ aVvD,t 


since the epitaxial layer is assumed to be uniformly doped with a density 


N In this model we have neglected the following physical processes: 


bh 
(l) “The Hite of the impurities due to electric field created 


during the diffusion process; 


(2) The growth "tail" and out-diffusion of donor impurities 
from the substrate that may occur into the epitaxial layer during the 


epitaxial growth; 
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(3) The precipitation of the impurities that may occur during 


the diffusion process; 


(4) The interaction of impurities with each other. 


2.4 Physical Description of the Diodes Selected 


Two basically different structures have been selected to 
represent the two distinct kinds of varactors. The choice of different 
structural parameters is based upon the discussion in Section 2.2. The 
doping profile is predicted by the mathematical model given in section 
2.3. The diode chosen to represent varactor diodes is named ‘diode No.1', 
and the diode chosen to represent step-recovery diodes is called 'diode 
Nos 2‘. Figures (2.1) and (2.2) display the doping profiles chosen for 
diode No. 1 and diode No. 2, respectively. Table (2.1) gives the 
description of the various parameters selected for our mathematical model 
to predict the doping profiles displayed in Figs. (21) “and: (@77 es tape 


(2.2) gives the description of the structural parameters of the diodes. 
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NET IMPURITY CONCENTRATION 
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Fig. 2.1. Net impurity concentration of diode No. 1 versus distance. 
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Fig. 2.2. Net impurity concentration of diode No. 2 versus distance. 


Ly 


°¢ “ON HGOIC ANV T °ON HGOIAa 
Od GHLORTHS SNOISNAWIG TVOISAHA °7°Z ATAVL 


SUOIOTW SUOIOTW 
g 0S ¢ 
SUOLOTW SUOIOTW 
01 0z 


ssouyoruL 
azokeT Terxeqrdy 


“woe 207 X 61°0 


"wD e 
Soe Oe 27610 


ssouyotulL 
ajerqysqns 


eso JO ePerY UOTIDeS-—ssoI9 H 


°"¢ “ON HGOIC GNV T °ON HCOIC AO Ssaooud 
NOILVOITSdVa AHL YOL GCHLOATAS SUMLENVEVd “1° AIAVL 


See 
uor0g snozoudsoug Ee 


TerzsjzeW 
1oqda.0y 


B8utdoqg 
TekeyT 
Terxeqyridq 


OUTIL 
uoTSNIFTG 


eianjerzodusy, 
uoOTSNFIITC 


TeTz9IeyW 
z1ouog 


eer) ie : - = 
a a 7 ; Lo » 
ae [tax 
a Mi = } 7 
wv) 
I 
, 
Bs 


Apes Lece 


A" 


q 


; etal 
+2 pk 


ise 


> EF SLED 


ar? 
> hh 


Jol 


— 
ai 


% 


+ Seve 


Ba 
“ey 


18 


CHAPTER III 
PHYSICAL MODELLING 


In this chapter, important physical processes occuring in the 
device are described, Various simplifications and assumptions are 
Stated. Mathematical equations modelling the simplified physical 
processes are formulated. At the end of this chapter, the validity of 


the one-dimensional model is discussed. 


3.1 Introduction 

As stated earlier in Chapter I, numerical techniques have been 
preferred over analytical techniques to avoid most of the inherent 
limitations in the latter approach. With the aid of numerical techniques 
it is possible to include the various physical processes in their most 
general form. However, many serious difficulties are encountered in 
formulating a general numerical scheme. These difficulties arise from 
severe non-linearities inherent in the mathematical equations describing 
the problem. The other limiting factors which are to be taken into 
consideration are the availability of a limited computer memory and the 
need to keep the computation time within practical limits. In view of 
the above considerations, various assumptions have been made to simplify 
‘the problem mathematically. In the following sections, these simplifying 
assumptions are outlined and a qualitative discussion is offered on their 


validity in the present case, 


‘an ! ft \eye SR toe = "a ah a 7] “ a4 Livy.” 
ma y Da. a al) 1 cA 7 i rt i 

os : a r) it i : : on 
ae - a 

: 7 Le j ; ¥ 1 
y i. i 
v oe | a = ’ 
; ea ae ‘ 
i : 
a) 

a 


4 ; TT) MAP OD ; 


POAT ee ih T Pia 


’ , ie 
ashy wi 7) eee . Meer os i 2. I rege! ee Gees nde a 


Sto &HBIIiq@ERey bor <4 DAS rah etrad sal tadhonel 
j 


hones Paste, 


Rajeyig BIAS Loe NeT rela 


0 


70 CeiLb Lis O15 \.SaToetio 91.0) "©: nee ee is me BAe 8 not ae 


To8itb @2° loner d jjnte 


ned. Sia Geary heise ei r soUnku ul ee | Laue whale § a 


e 

No te ty 9m tear biove 0 eaunreneieh Dp ii hee si. 
eotigraniee) Inbinonee 16 btw ately | ésbongymis n ate nt tind 
van a | 280K ten? pet eoanuwiny Teo apyalg wheldgny ‘cod Shu lone oe ate 
an Pi | ie bavrasaironie pon eabtbret2ith eqePybe ion , teva ured Pe 


ae Bra Serving? t1i5. sedi rede Ls Sitomue ty (one ~~ 


ee aiub bitin t fasttineotiu ods 4) whos ail. elke iw~ 


- oe eat nace en oY te Nolte fiocvont ep otnkh yortia ie 
ot en gen J ehikiidn boy tel pe Ao° yuk Ye Cheb a ee 
i sbi, Ae ciel RES ODA piiw agit aotan THUGS erke —_ bd. 


i) Libtite: ih i i Arad SRT 2 rae 25 Aha rneiiatact ail caged ‘aevede. Mt 
Are yh: i" 


} Mare a ; 
ae: 


Oe Pes 


My uy means s te roi yiitip aot atten Died ot ot): atria mere 


i 
oe tie, tid 
ry 


Us 


on | no beorsts bs dpaaepeit, Byles tay @ fs eb, Shield Stped ie OS ces 


RP an te | | es nee cits, wt ead 


ey i Ae geo A oie 
ee Oty Oy ee tas, 


19 
3.2 Important Physical Phenomena 


In this section, various assumptions made to simplify the 
mathematical problem are outlined, Mathematical equations to bepresent 
the simplified physical processes are formulated. 

Dezel bulk Generation ~ Recombination 

The measured properties of silicon p-n junctions at room 
temperature deviate considerably from that predicted by enoeiiey ae 
yideal" theory. Later, Shah, Noyce and Srockie, | improved the theory by 
taking into account the appar ey eae and generation of the carriers in 
the space-charge region, They used a model of the single-level uniformly 
distributed Shockley, Renae and Ene recombination-generation centers 


to-explain the p-n junction characteristics. According to Shockley and 


Rend °: 
pn - ey 
U SSR rR a i ga eee aoe aoe) 
(nak ay ss a (ped PLT, 
where 


. : -1 -3 
U = net recombination-generation rate (sec cm ~) 


-3 
p = density of holes (cm ~) 
| Eo 
n = density of electrons (cm ~) 
, Z -3 
n. = intrinsic concentration of carriers (cm ~) 
ae 


= life time of excess holes in heavily doped n-type 


po 
material (sec) 
ag = life time of excess electrons in heavily doped p-type 
no : 
material (sec) 
Pie density of holes when the Fermi-level is precisely at the 
| 3 
trap level E. (cm ~) 
nj = density of electrons when the Fermi-level is precisely at 


-3 
the trap level E, (cma) 
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Recently Deleon 7 calculated the effects of the spatial 


dependence of recombination centers on the I-V characteristics of p-n 
junctions. He showed that these Borers may explain the experimental 
deviations from Shockley ts! theory and from the theory of the carrier 
generation-recombination in the space-charge region. The 'gradual 
capture effect' is considered in his analysis by noting the possibility 
that the concentration of the capture centers tends to decrease when 
going inwards from the semiconductor surface and, therefore, for the 
Carrier lite time to increase in this direction. 

Equation (3.1) is derived under the assumption that majority- 
carrier density is large compared to trap density so that we may neglect 
the change in charge density produced by changing concentrations in the 
traps. Two other assumptions made are non-degenerate conditions and 
that the readjustment time for a trapped electron is negligible compared 
to the average time for the trap to emit the electron and capture a hole. 
In order to use Eq. (3.1), we have to know trap level EY - E and carrier 
life times. alguna obtained them for his theoretical calculations 
by matching the theoretical formula, using linear-graded junction 
approximations, and the experimental data at three points. In this 
approach, however, no such elaborate effort is made. Rather, these 
parameters are chosen from empirical data in the literature. As the p-n 
junction is formed by diffusion of shallow impurities, phosphorous and 
boron, EY - E can be taken to be 1okr??, The life times of holes and 
electrons, a6 and Ue respectively, are assumed to be equal and equal to 


07° sec, unless otherwise specified. A higher value, though not 
O eeierice: is used to decrease the effect of recombination which 
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should be the case in the ideal step-recovery diode. Moreover, we 
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treat carrier life-times independently of carrier densities, as the 


effect of such variation will be negligible in most beige * 


6.2,/ Mobiiity 


The dependence of mobility on doping density and field strength 


are well ete eae 


. §or simplicity, the carrier mobilitres ane often 
assumed invariant within a particular sample of semiconductor material. 
This, of course, is sometimes an unreasonable simplification. Experiments 
show that within a given sample of semiconductor material, the mobilities 
93 : 
of holes and electrons undergo large changes ~. Such changes arise at 
large values of electric field and at impurity atom concentrations. For 
this reason, in this work, we consider the influence of both the electric 
field and impurity atom density upon mobility properties of holes and 
electrons. 
; ; : 36 

In our analysis we have used analytic expressions for hole 

and electron mobility dependence on doping density and field strength. 


Caughey et ave derived these expressions empirically from published 


experimental data. Doping dependence of mobility is of the following 


form: 
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‘ 
REQUIRED IN CAUGHEYS EMPIRICAL 


EXPRESSIONS FOR MOBILITY. 
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The field dependence of mobility is given by the following expression: 


, wee: ee Go) 


where 


i hoae low field mobility 


= Critical field -streneth 


CELE 
= 20,000 V/cm for holes 
= 8,000 V/cm for electrons 
|E| =tapso lute yalue-oL thetelecericer ield, 


Another effect which is important is carrier-carrier collisions, which 

further decreases the mobility. This effect is only important at high 

injection levels. The possibility of a reduction in carrier mobility by 

carrier-carrier scattering at high concentrations has been discussed by 
‘ 5 37 emo 

Fletcher , Hall”~° and Davies~’. They have found reasonable agreement 


between experimentally measured values obtained on germanium and the 


values obtained by the formula 
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Pera” 


ead ed eetric’ constant. (11.7 for) silicon) 


m1 = effective mass of electron (1.1m, forswsiitcom Gr 20%) 
a = effective mass of holes (0.57m, for silicon at 20° ) 


h = Plank'ts constant 


n = p = carrier density. 


ae chevs showed that hole-electron scattering is important at densities 
greater than Tepe eninge Since for the dynamic conditions to which the 
diode is subjected, carrier concentration of this order will be 
restricted to a narrow region in the vicinity of the junction, we have 
neglected this effect to reduce the computation time, 
2.245 Degeneracy 

In this study, Boltzmann statistics are assumed to be valid. 
This assumption simplifies the mathematical problem in many ways; charge 
Resort phenomena can be derived from Boltzmann's transport equation, 
validity of Einstein's relation can be assumed and Boltzmann's relation 
can be used. 

As has been discussed by eiocialay age the difference between 
Fermi and Boltzmann statistics is negligible when the predicted carrier 
concentrations are small compared to the effective numbers of states in 
the conduction and valence bands. The effective density of states in 
the conduction and valence bands is given by 


is 3/2 3/2 
eft T 
300K, ees 7G) 


N = 2,5 x 10!9 


povat T =300K, using Mogg = 0.49m for holes and Mogg = 0.98m, for 
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2.4 X 10!9 om-3 
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It is evident that degeneracy will first become apparent in the valence 
band. As discussed in detail by shockieye 5 the Boltzmann statistics 


are valid up to 0.1 N, i.e., hole density of 0.35 X Weenie and 


Te Bate for silicon. From Shockley's 


electron density of 0.24 X 10 
discussion, the error in Boltzmann statistics is less than a factor of 
two for predicted carrier densities less than about 4N; a hole density of 
3.8 X 101? au and electron density of 9,6 X Ke ime for) silicon. 
itsepoint is further illustrated by Fig. (3.1). It shows that classical 
approximation always overestimates the carrier density in the valence and 
conduction bands. 

Duos + out ace Recombination and the Surface Boundary Condition 

Surface recombination is characterized by a parameter termed 

surface recombination velocity, just as the bulk recombination process is 
characterized by the life time of minority carriers. It defines the 
average velocity of carriers which maintain the recombination process by 
their flow to the surface. This parameter is dependent upon the carrier 
concentration and is constant only in the case of small deviations from 
the equilibrium concentration. The velocity S depends essentially on the 
treatment of the surfaces. Only by careful etching may a low value of 
surface recombination velocity be obtained. If it is high, injected 
carriers will be lost to the surface, and leakage current will be 
generated at the surface, the direction depending upon whether the 


operation of the device exposes the surface to a condition of minority 


carrier flooding or depletion, A value of 100 cm/sec is, in the opinion 
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of the author, a reasonable value for carefully etched devices. It can 


be ehowal’ chat 


iC a a0 
2) or oat et Con /) 
O 
where 
S = surface recombination velocity 
C = the mean thermal speed 


Ro = the probability that a single carrier in single 
collision with the surface will be sent back to the 


bulk rather than absorbed by recombination. 


The boundary condition for minority carriersis given Pea (for a n-type 


specimen) 


us = qS (p - PY wel (oe) 


where p and Pp, are perturbed and equilibrium values of hole distribution, 
respectively. It is evident from Eq. (3.7) that the maximum value of 
velocity of surface recombination is C/2. 

In our analysis we have assumed that thermal equilibrium 
prevails at! thelicontacts, 4 .le.—§ carrier, density at the contacts is 
constant and equal to its equilibrium value at all voltages of interest, 
Physically this means that terminal current on the p-side is totally 
carried by holes and terminal current on the n-side is totally carried by 
electrons. At low and medium forward-bias conditions, the minority 
carrier concentrations reaching the contacts (as given by Eq. (3.8)) are 
very small compared to the majority carrier concentrations; so it makes 


little difference in predicting the internal carrier distribution and the 
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28 
terminal current if we use this assumption, which greatly simplifies the 


solution. However, at large forward-biases, the assumption that thermal 
equilibrium preyails at the contacts has to be modified according to 

Eq. (3.7), particularly on the p-side which is heavily doped and shorter 
in width compared to the n-side. 

Along with the condition of thermal equilibrium, it is also 
assumed that charge neutrality prevails at the contacts. These two 
assumptions physically mean that the contacts are ohmic. Again the 
condition of neutrality will be violated at large forward~biases, 

The assumption of thermal equilibrium at the contacts is 
augmented in step~recovery diodes by the fact that the doping profile is 
so chosen that the built-in electric field tends to localize the injected 


carriers near the junction and prevents them from reaching the contacts, 


3.3 Assumptions for other Physical Phenomena 


(1) Temperature is assumed to be constant at 300°K, 

(2) Net doping density is independent of time. 

(3) Impurities are completely ionised. 

(4) Effect of radiation is negligible, 

(5) Tunnelling is negligible, 

(6) Skin effect is insignificant. 

(7) Frequency of the a.c. signal poneddered is low enough to 


allow the semiconductor to be collision dominated, 


ea.4 Validity of One-Dimensional Analysis 


The fact that our analysis is one-dimensional should be looked 
into a little more deeply, to obtain some qualitative estimate of the 


effect it will have on our results. Many of the microwave varactors and 
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step-recovery diodes have a mesa structure, as Bioen Dn apie 0352). Ate 
first glance, the validity of a one-dimensional analysis appears to be 
doubtful for both the static and dynamic cases. Let us start by 
considering the case when forward-bias is applied to the diode and 
steady-state conditions are allowed to prevail. The hole current will 
flow down through the substrate to the ohmic contact, Fig. (3.3) shows 
that the stationary forward current has a radial component. Its lateral 
spreading is, however, restricted by the fact that, at a constant voltage 
the current density rapidly decreases if the effective base width becomes 
too large compared to the ambipolar diffusion length. It is because of 
the bulk recombination and the longer the path that a carrier takes, the 
greater is the chance that it will recombine. We also see that the 
greater the recombination, the less will be the lateral spreading, At 
high forward biases, carrier mutual repulsion and non-uniformity of 
temperature will contribute further to lateral displacement, but again 
the bulk-recombination will tend to reduce this additional radial 
component by heavier recombination and hence reducing the current density 
and number of carriers that are laterally spread. 

Moreover, in step-recovery diodes, the built-in electric field 
tends to localize the injected minority carriers in the vicinity of the 
intrinsic layer and thus tends to prevent the carriers from going deep 
into the substrate and from spreading there laterally. 

Now let us go one step further and consider the situation where 
a diode, initially in forward-biased state, is suddenly reverse-biased. 
Conditions in the substrate are different particularly in the vicinity of 
the p-layer. Since initially charge neutrality is absent here, a 


stronger spreading can occur by mutual repulsion of holes, Besides, the 
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Fig. 3.2. Enlarged cross section of a doubly diffused mesa varactor (not 


to scale). 


CARRIER 
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Fig. 3.3. Carrier zone in the diode substrate under steady-state, 


forward-bias operating conditions, 


Fie, 3.4. One-dimensional p-n junction model, 
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current is flowing as a drift current and with relatively low 
recombination compared to the steady-state case; it can make good use of 
lateral Serene Lateral spreading thus reduces the current density 

and a smaller voltage drop may result when compared’ to the one-dimensional 
model, 

When a diode is turned-on by a high current pulse the displace- 
ment current will initially be dominant and then gradually particle 
current will take over and mobile carriers will be injected across the 
junction; in this case, recombination will be quite significant. Thus, 
the lateral displacement of the injected carriers will be checked in the 
later part of the transient as in the steady-state case, 

However, the assumption of constant uniform temperature 
throughout the device when it is being subjected to high power operation 
is not valid. The temperature towards the center will be higher than 
that near the surface, because of imperfect cooling of the heat 
generated. Moreover, the surface recombination velocity has a finite 
value. This accumulation of charge at the Sota: when the diode is 
driven into high forward-bias condition will contribute to further 
degrade the one-dimensional model. 

Thus, from the above qualitative arguments, it is obvious that 
the values of dynamic characteristics predicted by our model will be the 
least that may be expected, The performance of the device will be 
degraded by various effects (lateral spreading, non-uniform temperature 
distribution and finite recombination velocity at the surface, etc, ) 

: ZOIse. 
which are ignored in our present analysis. Moll et al. give an 
empirical estimate of the allowance to be made to the recovery time in 


pt -i- on diodes predicted by the uniform one-dimensional model. He 
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Suggests that 10 picoseconds per micron of i-layer width adequately 
accounts for the discrepancy caused by various effects which the one- 
dimensional model ignores. It is obvious that for transition times of 
0.5 ns or more this correction is relatively small (8% or less in our 


case). 


3.5 Synopsis 

The Shockley - Read - Hall single-level model is taken to 
characterize the carrier generation-recombination in our device. Carrier 
mobility dependence upon doping concentrations and electric field is 
considered. The empirical expressions given by Caughey et al.>° are 
used. The validity of Boltzmann statistics is assumed. It is further 
assumed that thermal equilibrium and charge-neutrality prevail at the 
contacts. Constant temperature, time independent net doping concentra- 
tions, complete ionisation of impurities and no tunnelling are also 
assumed. Skin effect and radiation effects are neglected. Moreover, the 
semi-conductor is assumed to be collision-dominated. It is pointed out 
that the one-dimensional model fails to take into consideration various 


effects such as lateral spreading of carriers, non-uniform distribution 


of temperature etc. 
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CHAPTER IV 
MATHEMATICAL MODELLING 


In this chapter, the mathematical equations based upon the 
physical modelling discussed in Chapter III are set up. These equations 
are then normalised to avoid recurring factors, They are then simplified 
for both the time-independent (steady-state case) and time-dependent 
cases for the one-dimensional model. At the end of this chapter, the 


boundary conditions for these normalised equations are formulated. 


4,1 Introduction 

Having made the simplifying assumptions for the physical 
processes occurring inside the device, the next step is to outline 
various mathematical equations which are taken to govern the physical 
behaviour of the device. The system of equations used in this analysis 

| Laan Wi see 
was previously outlined by Van Roosbroeck ~ in connection with 
investigations of the flow of holes and electrons in semiconductor 
materials. Appropriate boundary conditions, necessary for the unique 
solution of the above mentioned equations, are easily obtained from the 
simplified physical model of Chapter III. This system of equations 
along with the appropriate boundary conditions form our mathematical 
model, which completely defines the mathematical problem. The solution 
of this mathematical problem will yield the behaviour of the device as 


originally desired. 
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4,2 Fundamental Mathematical Equations 
4.2.1 Charge Transport Equation 

The Boltzmann transport equation (our assumption of non- 
degenerate semiconductor makes the Boltzmann statistics valid) with 
various conventional approximations (Moll, pp. 62 - say and 
supplemented by quantum mechanical restrictions (Pauli-exclusion 
principle) is taken to describe the transport phenomena in our model. It 
yields electron and hole currents as a sum of a diffusion component and 
a drift component representing Ohm's law. Such a simple form of current 
flow equation yields the following equation for holes and electrons, 


respectively: 


Heo) = Te DUP, VECese) 


uh (ef) ep te, b) Cel) 


Jt) = qu @,t) ae; t)-Vi(r,t) 
+qD_(F,t) Vn(e,t) CET, 
4,2.2 Charge Continuity Equation 
The continuity equations for holes and electrons are statements 


of conservation of charge and areas follows: 
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where 


UC@,t) = net rate of change of injected charge carriers due 
> 


to recombination and spontaneous generation 
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and 


By (s,) = rate of generation of holes (electrons) in 


the volume due to photon or other exciting 


agency. (not considered in our model) 


4,2.3 Poisson's Equation 
It relates the electrostatic potential to the net spatial 
electric charge in terms of mobile carrier densities and the fixed net 
impurity atom concentration N, which is the difference between the donor 


and acceptor distributions Na and Na respectively. 


q 


v2v(z,t) = 4 [ace.e) pb) n@)| Pou ees) 


€ 
where 


N(r) = Nor) - N,@) 


4.2.4 Total Current Continuity Equation 
= 9 
J, @,t) = J 3) iF US en = es VY, t) toe, 6) 


VeS (r,t) = 0 ote 5 (4 08) 


Equation 4.6 expresses the total current density as the sum of 


the electron and hole particle currents us and J_, respectively, plus a 


displacement current term. 


4.2.5 Auxiliary Relations 
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where 
Do and a = electron and hole diffusion constants 


respectively 


Hy and te = electron and hole mobilities respectively. 


4.3 Normalised One-Dimensional Equations for Steady-State Conditions 

It is very convenient at this stage to express the relevant 
quantities in normalised form; the set of normalization constants is 
chosen with the criteria of achieving the highest simplification in the 
relations of interest, The lise of normalisation factors used is given 
Emotable (451). 

The fundamental equations may be written in normalised form for 


the steady~state case as follows: 
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= G12 
3, = 5,@) + 5,0 (4,12) 


with the auxiliary relations: 


d¥ (x) 


E(x) = - dx 
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N(x) = N {Od - Ni @& oieses (Akl) 
ey ne ee) | Res 
p(x) =e [>p(x) - YG] Bede cole) 
o(x) = ux) n(x) + p(x) Up Ox) ease (gL GD) 


where 1? and Ph are electron and hole fermi levels originally introduced 


by Shockley, 


4.4 Normalised One-Dimensional Equations for the Time Dependent Case 
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4.5 Boundary Conditions 


4.5.1 Steady-State Case 


Assumption of carrier equilibrium at the contacts yields the 


following equations; 
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~ (4.27) 
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Assumption of charge-neutrality at the contacts gives equations 


(4,28). 
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The most general relationships between the carrier and current 
boundary values may be specified; only their functional dependence will 


be indicated here: 
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For contacts of the ohmic type as assumed in our model, equation 


(4.29) assumes the following simple form: 
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oa (4200) 
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where 
a ae a 2 = equilibrium carrier densities at the 


nh 


external contacts, 


According to Shockley's definition of quasi-Fermi levels, we can 


write the following equation: 
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Condition of charge neutrality at the contacts can be rewritten 


in terms of quasi-Fermi levels and electrostatic potential as follows: 
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4.5.2 Time-Dependent Cases 


Four boundary conditions on the mobile carrier densities for our 
device may be specified. These may be characterised by the same 


generality of relations, as given in Eq. (4.29), which are rewritten as: 
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Mint) = ake (Oye 1,009 | 
pO, t). = | CORE, 1,09 | 
4 ~. (4.34) 
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For ohmic-type external contacts as assumed in our model, the 


boundary condition as given in Eq. (4.34) assumes the simple form: 


n(0,t) = a, 
p(0,t) = P, 
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A dependent boundary condition on the slope of the electric 


field is furnished by Poisson's equation: 
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For ohmic contacts, it reduces to the following form: 
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As an aside, the electrostatic potential ¥(x,t) and terminal 


voltage vi) may be computed by subsidiary relations: 


% 
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4,6 Summary 

The simplified Boltzmann transport equation is assumed to 
describe charge transport phenomena in our model. The charge- 
continuity, charge-transport, total current continuity and Poisson's 
equations are taken to govern the device behaviour. Various physical 
considerations of Chapter III yield the necessary boundary conditions. 


All the equations are normalised and reduced to one-dimensional form. 
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CHAPTER V 


ANALYTICAL FORMULATION AND NUMERICAL TECHNIQUES 


FOR SOLUTION OF STEADY-STATE CASE 


In this chapter the details of the numerical solution of the 
problem for the steady-state case are discussed. The mathematical 
simplifications employed to derive various equations are outlined. The 
iterative procedure for HOA nn is then given. Various numerical 
techniques used are also discussed. A procedure to estimate the 
consistency of the numerical solution is also given. At the end, the 
numerical algorithm is extended to compute the incremental capacitance 


and the series resistance of the diodes. 


Det Licroduct ion 

Steady-state solutions are required in this thesis for two 
reasons: 

(1) To study the static behaviour of the device; 

(2) To furnish the initial values for the solutions of 

dynamic cases. 

A self-consistent scheme for the solution of the junction device 
problem was presented first by Gamele who applied it to the steady- 
state solutions of junction transistors. This scheme has an inherent 
limitation (due to present limitations of available digital computers) 
in that it cannot be applied to the case where a large reverse-bias is 
applied. DeMari‘? exposed some of the limitations of Gummel's scheme and 
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recombination model, 

In this thesis, steady-state results are derived for a 
relatively general model, as discussed in Chapter III, using basically 
the algorithm furnished by Comet ae The algorithm is modified to make it 
more efficient at high forward-biases where Gummel's original algorithm 
tended to fail. (DeMact did not have this problem as his model did not 
include the recombination-generation term, though he claimed that the 
extension is straight-forward). At very high recombination-generation, 
however, the present algorithm also failed to converge. Expressions for 
the electron quasi-Fermi level, hole quasi-Fermi level and electrostatic 
potential, which are treated as independent variables (as in Devartts- 
derivation) rather than obtaining them through a set of identities (as 


done by Gummer), are used. 


5.2 Derivation of Reduced Set of Equations 
Equations (4.7) and (4.8) can be written with the help of 
Boltzmann's equations (4.14) and (4.15) as 
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The constants of integration os and ea are obtained by using 
the boundary condition of oe) and 2&9) at x = L. Substituting these 


values of C_ and C , we get: 
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Equations (5.9) and (5.10) give the values of the hole quasi- 
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Fermi level and electron quasi-Fermi level respectively, in terms of 
electrostatic potential, carrier mobility, and generation-recombination 
rate.. The electrostatic potential is furnished by the solution of 
Poisson's equation (4.11). The carrier mobility and generation- 

Pee eniimerion rate, which are a function of the above three parameters 


and are treated as independent parameters, are given by Eqs. (3.3), (3.4) 
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and (3.2). An iteration procedure as outlined in the following section 


is employed to resolve the various non-linearities of the system. 


9.3 Iterative Procedure of Solution 

An iterative procedure may be employed to cope with several 
non-linearities of the formulated problem, namely Eqs. (5.9) and (5.10) 
for ee and DWE) in terms of ¥(x) which is given by Poisson's equation 
(eae <4, 11)) “and ae Tne :9) and U(x); the last three quantities depend 
upon the first three quantities, which are taken to be independent 
variables. For purposes of setting up this iterative procedure, it is 
convenient to rearrange Poisson's equation in a more appropriate form in 
terms of a quantity whose value at the boundary is zero. The correction 
5 {J oi eas between the electrostatic potential distribution obtained by 
consecutive iterations (j + 1) and (j) is introduced (as done by cumin cs 
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if the terms of order (x) and higher are neglected. 


Equation (5.16) represents a seeond~order linear differential 


equation in the unknown 5 eee. 


(x), with boundary conditions of Eq. 
(5.14) if the remaining quantities are available after the completion of 
the eee iteration. 

The complete iteration scheme is shown in Fig. (5.1). The 


applied voltage V, is specified, a trial potential distribution is 


A 
chosen, and the absence of recombination is assumed to start the first 
eycle of the main iteration loop (labelled "j"). The electric field E(x) 
and carrier mobilities u (x) and uw (x) are computed as a function of 

E e ® (x) d (x) 
position. From Eqs. (5.9) and (5.10), the values of e P and e 
are then calculated. Knowing these values, Eq. (5.16) is solved to yield 
corrections 6(x) for electrostatic potential. After obtaining the 
corrected value of the electrostatic potential, the Gy cycle is 
repeated until the correction 6(x) becomes sufficiently small. Then we 
include the computation of generation-recombination rate U(x). This 
change in algorithm proves to be useful for the computation pores 
forward-bias conditions where in the presence of a large generation-~ 


recombination term, any significant discrepancy in the values of the 


starting values of electrostatic potential from the true values tends to 
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DIODE DIMENSIONS AND DOPING PROFILE SPECIFIED 


APPLIED VOLTAGE IS SPECIFIED | 
TRIAL SOLUTION FOR ELECTROSTATIC POTENTIAL w (x) 
IS COMPUTED 


| 


RECOMBINATION - GENERATION RATE U(x) IS SET TO ZERO. 


ELECTRON AND HOLE MOBILITY jn (x), 2 p(x) DEPENDENCE 
ON DOPING AND ELECTRIC FIELD IS COMPUTED. 


VALUES FOR £.Q.F.L. AND H.Q.F.L. ARE COMPUTED FROM Eqs. 
(5.9) AND(5.10) IN TERMS OF PREVIOUSLY DEFINED QUANTITIES 


POISSON Eq. IS SOLVED TO YIELD ON IMPROVED SOLUTION 
FOR ELECTROSTATIC POTENTIAL 


HAS 
DESIRED 
CONSISTENCY 
IN w(x) BEEN 


YES 


YES 


CONSISTENCY 


COMPUTE RECOMBINATION — GENERATION ae 
RATE U(x) | 


COMPUTE ELECTRON AND HOLE DENSITY n(x), p(x), TOTAL 
RESISTANCE Rs, TOTAL CHARGE STORED Qp AND Qn. 


Fig. ed Creat iterative scheme for steady-state calculations of 


diode parameters. 
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be amplified, hence endangering the convergence, 

As is obvious, the success of this iteration scheme is dependent 
upon the choice of the trial solution of the electrostatic potential. iB 
the discrepancy between the values of the true solutions and chosen trial 
solution is small, the solution converges to its final value in a few 
iterations and thus consumes less computer time. For step junctions and 
asymmetrical linearly-graded junctions, first-order expressions are 
available which can be safely used to start the Soret rom But, LOE 
doping profiles which are neither step nor asymmetric linearly-graded, 
analytical theory is of no great help. Some general approximate 
expressions based upon the assumptions: 

(1) the device consists of completely neutral and completely 
depleted regions, 

(2) absence of generation~-recombination and 

(3) constant mobility 
are worked out and solved numerically. 

Overstratten and eta also recently published similar 
expressions to be used as starting values. However, under high forward 
bias conditions, these expressions did not help and a choice based on 
experience seemed to work out better. 

For very high generation-recombination terms, the starting value 
of U(x) = 0 failed to be a satisfactory value. The solution did not 
converge, but rather went into oscillation. Though several starting 
values were tried, no satisfactory answer was obtained. This feature, 
however, did not limit our analysis, as due to higher carrier life times 
Pea Contac did not hit the teabi of critical values ol generation- 


verombination term, But this feature is another limiting factor which 
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was omitted by ine ee and Demise 


5.4 Numerical Techniques 


5.4.1 Numerical Integration and Differentiation 

For the steady-state algorithm, suitable numerical techniques 
need to be chosen to approximate the numerical integrations and 
differentiations of interest. In the numerical execution of such 
operations, an error is usually introduced which often arises from 
truncating an infinite series at a eonfeniert point. The differentiation 
and integration formulae adopted here are based on the Lagrange 
interpolation scheme. The first impression here is eee higher-order 

interpolation schemes to reduce truncation errors. But, we have to 
remember that higher-order interpolation may give misleading results when 
used with functions that are not well behaved. Such functions have 
higher derivatives having large magnitudes which fluctuate quite rapidly. 
For such functions a linear interpolation between the known values 

would have provided a better approximation than a third-order 

Pot nomial ie. Unless the higher derivatives being approximated are well 
behaved and small in magnitude, it is wise to use linear interpolation. 
Various functions in our analysis vary very abruptly in the vicinity of 
metallurgical junctions. Their higher derivatives are not well-behaved 
and possess large magnitudes. 

When an effort was made to use mixed formulae for integration 
using different orders of polynomials, completely unrealistic results 
were obtained which showed an oscillatory pattern. So, as a general rule 
linear interpolation was prefered in the vicinity of the junction, so 
long as a proper density of points, into which diode length was divided, 


was kept at a 'safe' level (determined experimentally). 
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5.4.2 Solution of Poisson's Equation 
Poisson's equation ubéarranged in the form of Eq. (5.16) is to be 
solved in terms of the correction for the electrostatic potential 
distribution. To solve this second-order linear differential equation 


(Eq. 5.16) in the unknown 54 4 Bye) with boundary conditions 
so) = 0 and Me Gay =O : Pe (Ga) 


a straight forward technique of approximating the analytical derivative 
with a finite difference scheme, followed by the solution of the linear 
algebraic equations is adopted to solve Eq. (5.16). Rewriting 

Eq. (5.16) after discretization, we obtain the following system of linear 


algebraic equations. 
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A, is given by the following expression: 
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| 3, + sas 0 0 
BA H2 
a2 Sa ae 5 
H H2 H2 
0 oa = B, Sk i 
He ne 2 
0 0 gee ee le 
fe 
where { 
Vp eae ch ‘ Gye vay 
B. =e a +e P 


n = number of sections into which the diode is divided. 


H = spatial mesh size. 


To solve this system (Eq. 5.18), which contains a tri-diagonal matrix, a 
direct method using Gauss elimination and back-substitution is used, 
though indirect methods such as solution by interation could also be 
used. Though pivoting is ot a ee es enhance the accuracy of the 
solution of the system, the nature of our tri-diagonal matrix in which 
the main diagonal is large in magnitude compared to the other two 
diagonal terms, pivoting is not required. Thus the above technique 


yielded a fast solution of the system of equations. 


5.5 On the Accuracy of Numerical Solution 


One of the most important questions about any numerical solution 


is its accuracy. In our case, though an estimate of truncation errors 


can be made at each step, yet estimating the overall accuracy picture is 


‘ee Cy 0 () 7 ' 


9 f at? mae Hl 


kh ‘st shall . ify onl avetesen lo va a Py 

} p } ’ ICA, Tet bERes * at a 
- 5! ees bs ' ’ fie 1h Blissey i | vy A “a a aw ten ba 
i ‘ . a 


i Poi bry cated iele AA q+ Jnl fe ahead gniee f 4a 
We VY. 
roi - eet oof Binen ocotascsodu? vd DoLine 2 or Tp cass Dat seri ani 

‘ eats io: gop ose pA Ssosiide in Vebwanipaen: wi are sovke elas ay 


pains 


Hi My iy al: wh eteey il GMa tipe ht! 200 1G esuen oa eae wit 10, HRs 


in ‘ Lc 1 

a) ; vi / ee Tat Cat isi rat iaes 
Wier, lhe : i : * bel noe 
ALAM AD ed aad’ oh ad beragnes sbuyieess Ti ary at satis L 
ee, a 4 Dn A 


H wee ¥ : : wal 
Se aoptgilss)  dveta ene? bs per 260 al ja norte | maR 


irae a | Ns i Me 
RAG TE pes ravaye gas 4g colton tam 


; | Tee oe mols) Loe), Lie zn stot i. e304 A. HO Cit 
ha : ' os (BOG BRE sak cl iota: ; i ekcell 
a tye Mees ane hidien 
Phd SS Vie SRORH) erty bap. To "siognt shit at 
‘abit iweiey ia ange $e Ri: ‘ kwaiota ae cy _ 
: : ae a, Mie RP 
a Bene “) ny ae ra a ‘Ve ome aco “i “@ ais 
wa sain in Mt aye ioe mel a lL Late } jn a ") wh 


ey ve 


Te a ae NA ae iN Pea 
y Ha al f i “ , a iiel 


} t. 1 R HM _ 7 ent AG hone he ” ea ) ey Hi ; 
Q ; i een a ri (Ep Pa! v 7 ny ag iy mt uh ; i Te ; t oh | i i 


" are 7, ,. a) oe a Phy a 


54 


rather complex because of the various other errors resulting from 
different sources and mixing with each other in a rather complicated 

way. Instead of evaluating the net overall error, we used the approach 
of Peart: ie@., tO test|the consistencywor the Solution, )various 
analytical expressions, whose numerical counterparts display very 
different errors both in nature and magnitude (i.e. errors arising from 
the subtraction of large quantities of nearly equal magnitude, truncation 
error etc.) are employed to calculate the same quantity and the results 
are compared. To test this consistency the following set of relations 
are derived. 


(1) Electric Field 


E(x) = - BO. Juathal9) 
E(x) = pe") oy 3 G3) Wis n(x It at + E(0) F.C. 20) 


L 
{p(x") ~ n(x") + N(x") t dt + E(L) cerab 
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E(x) 


(2)Electric Currents 
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SO, we see that the final consistency of ¥, p, and n can be evaluated 
first by comparing results obtained for the electric field from Eqs. 
(5.19), (5.20) and (5.21). The discrepancies can further be exposed by 
comparing the results for hole and electric currents obtained by Eqs. 

©. 23)and (5.24) and (5.26) and (5.27) with equations’ (5.22)"and (9.25). 
respectively. (These equations are used to evaluate the distribution of 
current throughout the analysis). The last comparison also shows the 


overall effect of the various internal inconsistencies on the terminal 


properties. 


5.6 Calculation of Incremental Capacitance and Series Resistance 
5.6.1 Incremental Capacitance 


The numerical algorithm developed in Sections 5.1 - 5.5 is 
extended to compute the incremental capacitance and series resistance 
of the diode. If Q is the total charge (of holes or electrons) in the 


device, the incremental capacitance is defined at V = ves as 
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Writing Poisson's equation in integral form in terms of electric field at 


the contact, we get, 


forse a G9 air NGok dx = E(L) - E(0) 


When the above expression is differentiated w.r.t. voltage (V), we get, 


fe) O 
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tieor.h.s. is the contribution,of the.surface charge at’ the terminal, to 
the total incremental capacitance. This term is neglected here, as it is 


only significant at very high injection levels. Hence: 
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or, writing in terms of incremental quantities 
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It is seen that, to determine the value of the Fe ea 5 capacitance at 
any bias voltage, two solutions are to be obtained for two slightly 
different values of applied bias. 
5.6.2 Series Resistance 
The value of series resistance is generally deiner as the 
resistance offered to the flow of conduction current external to the 
Space-charge region. For.the case of epitaxial varactors, series 


resistance in the epitaxial layer is predominant and is given by 


x} L 
Bees a p (x) dx aoe oo) 
O | ay. 

Pp and p, are the p-type and n-type resistivities which are 
functions of carrier concentrations and their mobilities. xy and xX, are 
the boundaries of the space-charge region. As stated earlier, our 
physical model provides no well defined boundaries for transition region. 
Kennedy and Nerden” in their numerical investigations on p-n junctions 
provide a definition of space-charge region which is in substantial 
agreement with other investigations of the linearly-graded junction 
problem. They assumed that the space-charge region is terminated at 
locations where the electrostatic charge density is one half its maximum 


value on each side of the structure. This definition is used in present 


investigations to define space-charge region for the calculation of 


series resistance. 


5.7 Summary 


in this chapter, the analytical formulation of the equations 


M75) & (4 .12)); describing the behaviour of the varactor diode under 
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steady-state conditions (as discussed in the previous chapter) have been 
outlined. The electron quasi-Fermi level es hole quasi-Fermi level © 
and electrostatic potential ¥ have been chosen as independent parameters 
Cas in Ceneiee WOCK)<. (ine generation-recombination and mobility 
dependence on electric field, impurity concentration and carrier-carrier 
@ateractions has beén considered. Equations (5.9), (5.10) and (5.16) are 
the final equations solved according to the flow chart given in Fig. 
fori). sequations (5.9) and (5.10) are derived independently (as done by 
Penatio), Equation (5.16) is written in terms of correction in the 
trial solution of electrostatic potential YY. The overall algorithm was 
changed slightly to suit the cases where the recombination and 
regeneration term is quite dominant. But, for very high generation- 
recombination values, the overall scheme fails as the starting value of 
zero for the generation-recombination rate is found to be highly 
unsatistactory.. A criteria to test the accuracy of the final results is 
also given. At the end, the numerical algorithm is extended to compute 


the incremental capacitance and the series resistance of the diodes. 
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CHAPTER VI 


ANALYTICAL FORMULATION AND NUMERICAL TECHNIQUES 


FOR TIME-DEPENDENT CASES 


In this chapter, we will discuss details of the numerical 
solution for various time-dependent excitations applied to the terminals 
of the diode. Three cases are discussed: 

(1) the terminal current as a function of time is known; 

(2) the terminal voltage as a function of time is known; 

(3) an equation relating the terminal current to the terminal 

voltage is known, 
The mathematical manipulations employed to derive a reduced set of 
equations from the fundamental equations, Eqs. (4.17) to (4.22), are 
outlined. The boundary conditions needed to solve the reduced set of 
equations are derived. Then, the overall iterative procedure used to 
solve the reduced set of equations is outlined. The discretization 
scheme is selected from various available discretization schemes. After 
selecting the discretization scheme, details of the iterative procedure 
for a current-driven transient are discussed. Then algorithms for 
solution of the other two forms of external excitations are outlined. 
In the next section, various sources of numerical error and their effect 
on the accuracy of overall results are discussed. At the end of the 
chapter, various trade-offs between different errors, to decrease the 


overall error without increasing computation time, are also briefly 


discussed. 
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6.1 Introduction 
Results are required for time-dependent excitations applied 
across the diode terminals with and without some passive elements 
connected to the diode to study the internal and external properties of 
the diode under dynamic conditions. Various forms of excitations are 
necessary to study different phenomena. To study the injection of 
charge when a diode is turned on from an initial condition of thermal 
equilibrium, a large step of current is a good choice of external 
excitation. Similarly to study the details of the reverse transient, a 
reverse voltage step large enough to sufficiently reverse bias the diode 
is applied to the diode initially biased in the forward direction. To 
study the device behaviour under large signal operating conditions, a 
sinusoidal excitation is employed. To study the terminal and internal 
properties of the diode when driven by the above mentioned excitations, 
an accurate numerical solution of Eqs. (4.17) to (4.23) along with 
boundary conditions (4.35 and 4.36), which are taken to govern the 
behaviour of our device, is required. Deters has reported a numerical 
algorithm for the solution of these equations for a simpler model of 
zero-recombination and constant mobilities. He successfully applied it 
to the simple case of a highly asymmetrical step p-n junction. Another 
form of solution of this system of equations is furnished by Scharfetter 
and ean they applied this approach to the solution of the large 
eyes : A 
signal behaviour of IMPATT diodes. In our case, DeMari's algorithm is 
extended to our relatively general model which includes the generation- 
recombination rate, and the mobility dependence upon electric field, the 
net doping density and the mobile carrier concentration. This extended 


scheme is then successfully applied to our selected two types of diode 
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i x9 a 7 a 
Siiuctures, namely: p -.n and Pal Sa een ne eS riche ce 


Ot? -Pperivation of Reduced Set: of Equations 


The fundamental equations (4.17) to (4.23) with the auxiliary 
relations Eqs. (4.24) to (4.26), which are taken to govern the behaviour 
under time-dependent excitation conditions, may be rearranged to a form 
more appropriate for numerical solution. If the expressions (4.17) and 
(4.18) for electron and hole current densities are inserted into 
Eq. (4.22), anduse is made of the auxiliary relation of Eq. (4.24), the 


displacement term may be explicitly written as: 


CE ty _ 


- = uy) p(x,t) + u(x,t) n(x,t) | E(x,t) 
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Moreover, if the expressions (4.17) and (4.18) for hole and electron 
current densities are inserted into equations (4.19) and (4.20)* 


respectively, we get: 
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and 


5 Sie (x; b) 
on({x,t) _ u(x,t) o*n(x, t) +{ —P_____ -_ B(x, t) BCs) “ROG 0) 


ot 52 ox 
du (x,t) 7 
OE 
-}- 2 wee, t) + BORD yy Gey | aye) 
oe (OG Gig a) as D) 


end, with the help of Poisson's Eq. (4.21), 


2 Ou (x,t) 
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ox 
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2 3 u CRs) 
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The reduced set of Eqs. (6.1), (6.4) and (6.5) with appropriate 
boundary conditions and initial conditions and external excitation J(t), 
as discussed in the following sections, represents a complete formulation 
In this formulation, the electric field E, the electron 


of the problem. 


density n and the hole density p are chosen as the independent quantities 
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and represent the unknownsof the reduced set of equations. 


6.3 Boundary and Initial Conditions 


For a two-contact device, four boundary conditions on the 


mobile carrier densities may be specified. 
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A dependent boundary condition on the slope of electric field 


is given by Poisson's equation (4.21) 
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For external contacts of ohmic type, the boundary conditions of 


(6.6) assume the simple form: 
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where n> Py Bs and By are the equilibrium carrier densities at the 
external contacts. The dependent boundary conditions become then 


9E(x, t) 
ox 


= 0 eee 6. 
x=0 aed) 
and are a consequence of the charge neutrality condition of the ohmic 
contacts. 


Moreover, thevelectrostatic potential Y(xjt) and terminal 


voltage Vv 6e) may be computed by the subsidiary relations: 
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where ¥(0,t) may be taken as reference value for electrostatic potentials 


and Va is the diffusion potential defined by 


L 
V.. = fos dx thermal sHeasul Onalicy) 
O 


equilibrium 


The electric field and mobile carrier density distributions 


throughout the interior of the device at the initial time t=0 : 
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the time dependent solutions starts at steady-state condition, as in our 
case, the initial conditions are available from the steady-state 


solutions obtained with the method described in Chapter 5. 


6.4 Iterative Method of Solution 

Equations (6.1), (6.4) and (6.5) describe a system of three 
partial differential equations in terms of three unknowns; the electric 
field E, the electron density n and the hole density p. The equations 
are in two-dimensions namely time t and position x. Equations (6.4) and 
(6.5) are non-linear second order partial differential equations of the 


47, 


parabolic type. Finite difference methods are chosen to achieve 
the solution of these two equations. The first step towards a numerical 
solution requires a discretization of relevant quantities at a finite 
number of points in both time and space. Appropriate numerical 
techniques are then to be devised to determine the spatial distributions 
of unknowns that satisfy the set of equations for assigned boundary 
conditions at each instant of time. An alternate approach of 
simultaneously determining the unknowns in both time and spatial domains 
first for broad mesh size and then gradually refining the mesh size till 
desired accuracy is reached, could be used. Because of heavier memory 
requirements in this approach, it is not considered feasible, 

The discretization problem for the partial differential equation 
of the parabolic type is quite extensively discussed in the next section. 
In this section a numerical method of the solution is sketched 
schematically. For our model Eqs. (6.1), (674) ‘and °(6.5); boundary 
conditions (6.8) and (6.9), initial conditions (6.13) and external 


excitation J(t) give the mathematical formulation of the problem. The 
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66 
expressions for mobility dependence on electric field, doping density 
and carrier density and thé net rate of generation—-recombination, as 
discussed in Chapter III, are also needed. The highly non-linear 
character of the equations requires an iterative procedure at each 
instant of time to reach desired consistency. From the knowledge of the 
spatial destributions at tabs it is desired to determine the unknowns 
at the next instant of time t=t,. An educated "guess"for the electron 


1 


and hole density spatial distributions at ty is made and these are 
inserted in Eq. (6.1) to solve for preliminary electric field spatial 
distribution at ty which in turn may be inserted together with the 


' for the hole density distribution, in Eq. (6.5) to solve for the 


"suess' 
electron density spatial distribution at t=t)- These improved 
distributions may be inserted in Eq. (6.4) and the iteration cycle at the 
instant of time t, may be repeated untill the desired consistency is 


1 


meached, at th The distributionsat time ty are then taken as initial 
conditions and the whole procedure is repeated to determine the 
distributions at time to. The process terminates when the final stage is 


reached. The scheme is summarised in Fig. (6.1). 


6.5 Selection of Discretization Scheme 
In general, the continuous quantities appearing in the 
analytical formulations are discretized at a finite number of points in 


both position and time co-ordinates; the ensemble of these points 


constitutes a "grid" or "mesh". Finite difference schemes are then 
employed to approximate the analytical differentiations and reduce the 


problem to the solution of a system of finite difference equations, which 


represents the discretized formulation. 
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6.1. Iterative scheme for c 
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The basic problem here is the selection of a discretization 
scheme featuring both convergence and stability, together with a 
satisfactory degree of efficiency of computation. Convergence is defined 
as the discrepancy between the exact solution of the original analytical 
formulation and the exact solution of the system of difference equations. 
The theory of convergence deals with the conditions under which the 
discretization error becomes equal to zero in the limit of vanishing 
distance between adjacent position and time points. In addition, the 
actual numerical results obtained as a result of computations on the 
basis of the discretized formulation, differs from the exact solution of 
the system of difference equations by an amount equal to "numerical 


error", The problem of "stability" deals with conditions under which 


the numerical error is small throughout the entire range of the solution. 


47,48,49,53 furnishes 


The available theory of numerical analysis 
several discretization schemes, the convergence and stability of which 
is proven for simpler cases. The theory is only complete for the case 
of linear equations with constant coefficients; for linear equations 
with non-constant coefficients, theory is still far from satisfactory. 
For non-linear equations, the conditions of convergence and stability 
are only available for a few special solutions and discretization 
Berens se Hence the only guide line available to us is the solution of 
the simpler cases from which conclusions are extrapolated for more 
complex systems of non-linear differential equations such as ours. 
Various discretization schemes available to us are discussed in 
fee cate foe As we see there, in general two types of discretization 


schemes are available: explicit schemes and implicit schemes. Stability 


and convergence requirements for certain types of discretization schemes 
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ece Cependent upon restrictions on mesh size. For a linear différential 
equation with non-constant coefficients, the stability and convergence 
are proven in Ref. 47 (pp. 107) for the simplest explicit scheme. The 


Bestriction on mesh size is 


< ia fw Peet Owls) 
(Ax) 2 2a (x,t) 
where 
2 
a(x, t) = coefficient of the ease). term. 
ox 


When applied to our case of Eqs. (6.4) and (6.5), we get: 


er ed 5 Ce AWAD: 
(Ax) 2 2u (x,t) 
5 ya eee Niles EGTA 
(Ax) @ 2u_ (x, t) 


Since u_(x,t) is greater than PBL the greatest value of 
n 
u_(x,t) will determine the critical relationship between space and time 
n 


mesh sizes. For 


Ax = 10° 
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W 


At 


This is the least value necessary to ensure convergence and stability. 
The value of At furnished by Eq. (6.15) is highly undesirable because of 


the huge number of points required for the time co-ordinate, resulting in 


a huge amount of computer time. 
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If we look at the implicit discretization schemes, we see that 
for the simpler case, unconditional convergence and stability without 
any constraints on the mesh size are guaranteed. This is proved in 
Ret 47e(pp. 102). for.linear equations with constant coefficients. 
However, as is obvious a partial differential equation of the parabolic 
type, discretized with an implicit scheme, requires the simultaneous 
solution of difference equations (one for each spatial point) at each 
instant of time, not required for the discretization generated by 
explicit schemes. But the freedom available in the selection of mesh 
size in both position and time co-ordinates appears to enable us to 
choose a reasonable number of points on the time co-ordinate and thus 
enable us to save computing time rather than increasing it. This is 
later verified in actual numerical computation. Due to this salient 
feature of Boe implicit discretization scheme, they are chosen over 
explicit schemes. The generalized implicit discretization has been 
employed (6=1), though the Crank - Nicholson discretization schemes 
(6=1/2) and other implicit discretization schemes (0 < 8 < 1) could also 
be employed. The Crank - Nicholson scheme has been demonstrated to be 
very sensitive to any numerical errors (errors in initial solution, 
truncation and round-off errors) and its stability is endangered even 
though the numerical errors are not severe from practical considerations. 

After selecting the discretization scheme, the derivatives are 
approximated by finite-difference expressions derived from Lagrange's 
interpolation scheme. We have to be careful in choosing the finite 
difference expressions used to discretize the derivatives. If the order 
of the difference scheme is higher than original differential equation, 


the approximate solution may contain extraneous solutions which are not 
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at all related to the true solution. These extraneous solutions 
severely affect the convergence and stability conditions. The case of 
an ordinary linear differential equation with variable coefficients has 
been treated in reference 51 and the presence of extraneous solutions is 
found to be very hazardous, Although higher order interpolation schemes 
seemed to improve the accuracy of the numerical solution, their use is 
avoided in this work due to their highly uncertain stability and 


convergence conditions. 


6.6 Details of Adopted Procedure 


After choosing the generalised pure implicit scheme for 
descretization for the solution of Eqs. (6.4) and (6.5), we proceed to 
discuss the details of the iterative procedure outlined in Section 6.4. 
Bauations (6.1), (6.4) and (6.5), when discretized on the basis of the 


above scheme, generate the following system of difference equations. 
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where 
Ou 
(AGES aa 
k Ade fe ik ot 
Big = 
an Yd ea 
od rn 
jet vee 5 on 
At ~ kind Soe as ees Ox 
2 at 
k Be af 
a rs = 
i: ay (2% kPa ",) k nt , ie 
CRED) 
where 
OU 
nN 
Co Ox ete eh 
K i 
pera pale 


ine left hand index (k = 1; 2;°3,4... .\£) represents the time instant 
in the discretization scheme. The right hand index (i = 1,2,3,....2) 
represents the position in the discretization scheme. Quantities at 


instants of time t < and at any spatial point 0 < x <li, ae 


14 
Supposedly known and the solutions for the next instant of time ,t are 
sought. According to the method briefly outlined in Section 6.4, an 
iterative procedure is employed at each instant of time to determine the 
unknown quantities. Toystart the iteration at each instant of time, 


trial solutions have to be selected. The values of hole density, 


electron density and electric field computed at the previous instant are 
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found to be satisfactory to serve as the starting solution. 

It is convenient for us to introduce the corrections An and Ap 
for the electron and hole densities, respectively, defined as the differ- 
ence between densities corresponding to two successive iterations 
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where 


=O Oe (number of iterations) 


If the convergence of the approximate spatial distribution 
occurs during the iterative procedure at each instant of time, the 


following relations are readily verified with the aid of definition 


(6.19). 
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In addition, the time dependent property of the boundary conditions in 


fig. (6.8): yields: 
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If the quantities labelled k in equation (6.17) and (6.18) are 
identified with the subscript (j+1) (with the exception of the terms 
Pi and e124)? the relations in Eq. (6.19) becomes: 
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Moreover, the electric field at the instant (k+l), may be 


explicitly recovered from Eq. (6.16) and rewritten as: 
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If the second order terms (Ap)2 and (An)? are neglected and the 
three point formulae are used for numerical computation of derivatives 


@£ 4n and Ap in Eqs. (6.18) and (6.19), these assume the final form: 
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The numerical derivatives of the electron and hole densities in 
the expressions FP and FN are computed from the two-point differentiation 
formula, since, due to the steep gradient at the metallurgical junction 
of the relevant quantities, numerical derivative formulae based upon 
higher order interpolation formulae yield highly inaccurate values in 
the vicinity of the metallurgical junction and in some cases exhibited 
an oscillatory pattern for the distribution of internal quantities. The 
general rule followed in Chapter V; that for a function which does not 
exhibit well behaved higher order derivatives, lower order interpolation 
yielded better results, is observed. The only point about which we need 
to be very careful is that we choose a high enough density of points in 
the space co-ordinate so that truncation errors are bounded to ensure 
numerical stability. 

Equations (6.27) and (6.28) may, with the help of boundary 


conditions of Eq. (6.21), also be written in the vector form as follows: 
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ak (any DP 2 Ge Re 200) 
k 
where 
a and Th are tri-diagonal matrices 
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COMPUTE IMPROVED VALUE OF CARRIER MOBILITIES FROM 
EQ. (3.4) AND(3.5)(IT CONSIDERS MOBILITY DEFENCE ON 
ELECTRIC FIELD AND CARRIER DENSITY) 

¥ 


USE EQ. (6.27) AND n(x,t), E(x,t), Up(x,t),U(x,t) AT tp TO 
SOLVE FOR IMPROVED HOLE DENSITY p(x,t) SPATIAL 
DISTRIBUTION. 


Moe EQ. (6.28) AND p{x.t}), E 
SOLVE FOR IMPROVED ELECTRON DENSITY n(x,t) SPATIAL 
DISTRIBUTION. 


USE EQ.(3.2) AND plx,t),n(x,t) TO COMPUTE 
IMPROVED VALUE OF RECOMBINATION-REGERATION 
RATE Wi(x,t)¢ 
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Fig, 6.2. Continued. 
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As weysee, Eqs, (6.26), (6.2/7) and (6.28)3 are our final 
equations on which our whole iteration scheme is to be based. 

First the structure parameters: diode thickness, junction area 
and diping profile are specified. Carrier mobility dependence on net 
doping density and electric field is then incorporated (given by Eqs. 
(3,3) and (3.4)). Then the formula for calculating the rate of net 
Pe tio Mpbadihd is added (given by Eq. (3.2) in our case). The 
Steady-state distributions obtained in Chapter V are specified as the 


initial spatial distributions of electric field, and mobile carrier 
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meeCilFY DIODE STRUCTURE PARAMETERS. ( JUNCTION » AeA. | 
_ EERE SN ERS, PD Byte NG FROGILE) | Ut ae 


BRECIFY THE INITIAL. CONDITIONS AT t=t=t........... 
Esc Fo:) 
cr Osaka.) ALSO CU px sto) Unis eke ilto) 
ni xetos 
SPECIFY THE BOUNDARY CONDITIONS FOR MOBILE FOR MOBILE CARRIERS 
DOr )e= 1 ce Wa 


COMPUTE SSC OGT Tite er FQ.1 23) (WHICH Gives 
MOBILITY DEPENDENCE UPON NET DOPING DENSITY) 
SPECIFY THE EXTERNAL CURRENT EXCITATION | I(t) | zo rm 


oe : ee 
GO TO NEXT TIME INSTANT 
Pea tel Coup Ke lee Py 


SPATIAL DISTRIBUTION n(x,t), p(x,t) AT ty; ARE TAKEN 
TO START ITERATIVE PROCEDURE AT ty, 


USE EQ. (6.26) AND ate tp VAL TIME fet esl xet) ticy 
Up(x,t), AND Un(x,t) AT t, TO SOLVE FOR ELECTRIC FIELD 
E(x,t) SPATIAL DISTRIBUTION AT ty 
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fae. 6.2, Detailed overall iterative procedure used for the current- 


és driven transient case. 
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82 
densities at the instant k=0, An external excitation of total current 
as a function of time is specified to drive the evolution of the 
transient. The spatial distributions of E, n, and p, the unknown of the 
problem, are sought at the first time instant k=1 and will be obtained 
with the use of an iterative successive approximation scheme. The 
iteration index (j) is set to zero and the spatial distributions of the 
mobile carrier densities available at the previous instant k=O are taken 
as trial functions to start the iteration at the instant k=l. Equation 
(6.26) yields the spatial distribution of the electric field which may 
be inserted in the implicit system of Eqn. (6.27) to solve for the 
corrections soe, oe Ui =e2, S40 © . &-l) Tor hole density. mauacion 


(6.23) then yields an improved hole density distribution which, together 


with the electric field distribution last obtained from Eq. (6.26) 


(j+1) 
BAL ; 


improved electron distribution. Then, Eqs. (3.4) and (3.5) are used to 


solves for the corrections so that Eq. (6.22) yields an 

yield an improved value of mobility based upon electric field and carrier 
densities. Equation (3.2) is used to yield ot, daaeete value of the net 
generation-recombination rate. The terminal voltage is obtained by 
integration of the electric field. Then over-all consistency tests are 
performed, A muetite outcome of these tests at the considered instant 
of time (k=1) increases the iteration index by unity and restarts the 
cycle with the insertion of the improved mobile carrier distributions in 
Eq. (6.26) for an improved electric field distribution at the same time 
met ie (k=1). A positive outcome of the consistency test accepts the 
quantities obtained by last iteration as 'the' solution at the time 


instant considered. Unless the final state (KEE) Ws reached, the time 


index is increased by unity, the iteration index (j) is set to zero and 
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83 
the solution at the next time instant (k=2) is sought. Ihe iteration as 
restarted by taking the mobile carrier distributions and mobilities 
available from the previous instant (k=l) as the trial function and 
inserted into Eq. (6.26). The same iteration procedure is repeated to 
achieve 'the' olor ae the instant (k=2) and the subsequent instants 
Soucime until the final state (k=f£) is reached, 

The tri-diagonal matrix system is solved by the same techniques 
successfully employed in Chapter V for the solution of the tri-diagonal 
System generated by the discretization of Poisson's equation in the 
steady-state scheme and is described in reference 47 (pp. 103 - 105). 

The overall consistence, between the unknowns E, n and p and 
the system of equations determining them, may be tested at each iteration 
cycle for each instant of time by considering the magnitude of the 
Pectections An, and Ap and the difference of ,the electric field 
distributions obtained by two successive iterations. If the magnitudes 
of these quantities are smaller than any specified value (iteration 
error at each instant of time), achievement of overall consistency at 
that instant is assumed. The final state (k=f) is reached when the 
external excitation has reached the final state, and the unknowns E, n 
and p are unchanged for an increase of time. 

A most demanding requirement for the success of the described 
method is the convergence of the iteration procedure to the overall 
consistent sclutions from a set of approximate trial distributions at 
each instant of time. In all the cases tested for various structures 
(two more than reported here) and external excitations, convergence was 


always obtained, though under certain conditions stability was 
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endangered (weak stability). The rate of convergence was quite fast in 
situations in which significant variations (in time) of the total current 
occurs, so that the electric field (given by Eq. 6.26) becomes the 
quantity driving the variations of mobile carrier distributions (Eqs. 
6.27 and 6.28) and is essentially responsible for the evolution of the 
transient. A slower convergence rate was observed under conditions, 
opposite to those just described, i.e., when total terminal current 
remains constant. The rate of convergence for such cases was found to 
be well within practical limits except for conditions arising towards 
the termination of the transient, when the quantities approach the final 
steady-state distributions. The final distributions could be obtained 
from the steady-state schemes using current as the driving quantity. 
However, this was not the limitation in our case, as we aim to see the 
phenomena for comparitively shorter times. These effects were 
comparitively similar to that obtained by Denaciet for his“léess general 
model of constant mobility, and zero-recombination. But, the convergence 
was not that fast in the case of the ‘electric field driven' transient 
because of the dependence of the carrier mobilities upon the electric 
field, which therefore had to be adjusted to reach overall consistency. 

Tt is clear from the method and formulation described that 
complete freedom is ai ee in the selection of both the spatial and 
time steps. Their selection depends upon the total error that can be 
tolerated. The time-step distributions may be conveniently chosen and 
automatically adjusted by the computer during the evolution of time. 


This is discussed in Appendix A. 
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6.7 On Different Forms of External Excitations 
One of the following three possibilities is always known at the 
terminals of the diode: 
(1) The terminal current as a function of time (current—driven 
case)* 
(2) The terminal voltage as a function of time (voltage-driven 
case)3 
(3) An equation relating the terminal current to the terminal 
voltage, 
6.7.1 Procedure for Voltage-Driven Case 
So far we have discussed the current-driven case. Now if the 
applied voltage as a function of time is specified at the terminals of 
the device, we need to substitute Eq. (6.26), expressing the variation 
of the electric field in terms of total current, by Poisson's equation 
(Eq. 4.21), which is written in terms of electrostatic potential ¥(x,t). 
The applied voltage then appears explicitly as a boundary condition for 
Y in the solution of the second~order differential Eq. (4.21). Now Eqs. 
eee) and (O65), along with Eq. (4.21), rewritten in terms or Che 
variables ¥, n, and p, comprise a new set of fundamental equations. Eqs. 
Bor rand (6,5) are treated just as for the current-driven transient 
case, Equation (4.21), the new equation of the set, is discretized and 
solved by the same method as was used in the solution of Poisson's 
equation for the steady-state case. An iteration procedure of the same 
type as described for the current-driven transient is used to reach the 
desired consistency between ¥, n, p and Eqs. (6.4), (6.5) and 4521)" 3b 
each instant of time. A slower rate of convergence of the iterative 


procedure was observed in this case compared to that of the current- 
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driven transient case. This is expected, since the 'error' due to the 
variation of the electrostatic potential at the boundaries is contained 
at the boundary (at one contact) at the beginning of the iteration scheme 
at each instant of time and propagates throughout the whole interior 
before the solution is allowed to proceed to the following time instant. 
Moreover, the solution of this system of equations requires a heavier 
computation load, since it involves the solution of an additional tri- 
diagonal matrix necessitated by the solution of Poisson's equation. 

An alternate approach is to use the method suggested by ene 
for his simpler model. In this method, the current-driven transient 
solution i.e., specifying the total current J(t) and determining the 
terminal voltage vo at each instant of time, is used as an integral part 
of the solution for the voltage-driven transient. The solutions for 
total current and internal distributions are assumed to be known for 


eimes ty < and quantities at the following time instant are sought. 


ii. 


If the total current i were known, the iterative procedure described 


in Section 6.6 could be applied to obtain the unknown distributions at 


the instant wee The problem now consists of searching for a particular 


value of current Ko that will yield the specified terminal voltage ie 


This is achieved in two stages using the current-driven transient 


algorithm at the specified instant of time. A "nrediction" stage 


Gye a 
initially furnishes two estimates of the total current k J an a 


close to the 'correct' value 1? which is to be determined. This is 


achieved by using linear (two-point) and parabolic (three-point) 
extrapolations on the current responses as a function of time which 
Two independent 


furnishes, in general, two estimates for the current i 


applications of the current-driven transient algorithm yields for the 
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Fig. 6.3. Prediction and interpolation stages as proposed by DeMari 


for computation of the voltage-driven transient case, 
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(1) (2) 


two currents k J and ic J the corresponding values of terminal 


1 2 
voltages - y and = ve which, in general, differ from the specified 


terminal voltage iva: A second stage based on interpolation of the curve 


(h) Ciel oee nerves oo) as a function of ie 


3 a J at the same time k, 


may now be started to obtain an improved estimate of the total current 


that generates a terminal voltage within an assigned tolerance about 


the specified value. A linear interpolation on the graph of hoy 


(2) (1) 22) (3) 
2 ee ioe 


k 


J versus re and Be yields a third estimate 


corresponding to the specified terminal voltage ka’ A third application 


of the current-driven transient algorithm at the instant k for the 


(3) 2) 20 0K3) 


current o J furnishes a third point aoe 


voltage graph, so that a parabolic interpolation may be used to further 


ve on the current- 


improve, if necessary, the estimate of the total current. The procedure 


(h) 


continues for higher-order interpolation schemes until values k A 


(h) 
ae 


VY ,» is determined. Solutions at the instant k are then available and 
a 


the overall procedure may be restarted at the following instant ktl with 


that generates the voltage sufficiently close to the specified 


the prediction stage. For a few initial time points (k = 1, 2) the 
prediction stage may not be performed as described; an arbitrary 
estimation procedure may be employed equally well, if the time step is 
kept small enough, in order to obtain a quasi-linear current-voltage 
characteristic in the interpolation stage. For a large oleae step, 
this procedure tended to fail, due to the large variation of current; 
the interpolation and extrapolations may fail to give proper values, 
unless an extremely small time-step is employed. In this case the first 


procedure using the electrostatic potential, the hole density and 


the electron density as independent variables is recommended 
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as there is no chance of failing. However, for snald voltage driven 
transients, the 'extrapolation-interpolation' procedure is very 
efficient. The overall scheme, using this procedure is outlined in 
Proure (6.4). 
6.7.2 Outline of Procedure When an Equation Relating 

Terminal Current and Voltage is Known 

When a general two-port network, composed of ideal current and/ 
or voltage sources and passive circuit elements, is connected to the 
terminal of the diode (eons device), both the total current 
density J(t) of the device and terminal voltage Vf) are quantities 
dependent upon the properties of both the diode and the network. At a 
particular time instant k, the two unknowns ,J and Ve are determined 


k k 


by the system: 

J) oas(Oe 49) 
oF eG. 40) 

io 5 Ove 


The above equations are the discretized current-voltage 
relationships at the time instant k of the diode and the network 
respectively. These equations can be directly incorporated into the 
iterative procedures described previously. For a reasonable number 
of elements, the time required for the solution of this system is 
insignificant compared to the time required to compute one iterative 
cycle for the diode. The undesirable effect is the resulting slower 
convergence rate, so that more iterations are needed to achieve overall 


consistency at each instant of time. 


For the simpler case in which only a resistance R in series with 
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FIRST PHASE: PREDICTION ON THE 
J=J (t) CHARACTERISTIC 
(1) (2) 
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Fig. 6.4. Overall scheme for computing diode parameters for the voltage- 


driven transient case. 
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91 
a time-dependent voltage source VQ Ct) is present in the external Crreunts 
equation (6.46) results in: 
es, 
; NOs 
Njaddtwik! Gick’ aun Ne 


ASA 8 2. (6.47) 
. RA J 


where 


A = area of the device 


Ve = thermal voltage kT/q 
J = q _ eet. 
(@) SanEENEReEanEneaiemnemene eel 
LD 


and voltages and currents are in normalised form. Results have been 


computed for both step and sinusoidal voltage sources. 


6.8 On the Accuracy of the Numerical Solution for Time-Dependent Cases 


In the time-dependent solution discussed in the previous 
sections, truncation errors both in the spatial and time domains are to 
be considered. In addition, the interaction between the errors 
introduced at the different time instants becomes a dominant factor for 
Several. accuracy, The discretization error in the time domain and) the 
growth of numerical error for increasing time are responsible for this 
interaction and consequently determine the conditions of convergence and 
Stability, respectively. Due to this fact, the conditions of convergence 
and stability for different external excitations will be different and 
will have to be tested separately for each case. In the following sub- 
sections, the principal sources of error arising in the time dependent 
solution (besides the inaccuracies in the spatial domain), together with 


techniques considered suitable to evaluate and control the influence of 
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92 
such errors on the final results, are briefly discussed. 
6.8.1 Truncation Error in the Time Domain 

In the formulation adopted for either the current-driven or the 
voltage-driven transient, the time derivatives are approximated by a 
discretization scheme, resulting in discretization error. This 
approximation is often achieved by truncating a Taylor expansion at a 
convenient point to obtain simple finite difference expressions. This 
truncation error depends upon the values of the higher derivatives, on 
the size of the time step, and on the order of the finite-—difference 
equation used. The magnitude of the higher order time derivatives of 
the functions to be numerically differentiated is very much related to 
the abruptness of the time variation of the unknowns, directly dependent 
upon the time variation of the external excitation, for a given time 
step. For a decreasing time step, the variation of the electric field 
is more emphasized, if compared to the variation of the mobile-carrier 
densities. This is a consequence of the physical properties of the 
mobile carrier densities (and therefore the particle currents) 
requireing a finite time for a significant variation, where as the 
electric field (and therefore the displacement current) may respond 
immediately to any external excitation. 

Time-mesh size compression when the external excitation varies 
more abruptly during the transient (resulting in larger values of the 
higher derivatives) to keep the truncation error constant and within 
limits can be employed and is discussed in detail in Appendix A. 


One tempting idea is to employ higher order finite difference 


schemes to approximate the time derivative, to reduce the truncation 


errors. However, the stability and convergence may not be guaranteed 
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for each case. In some particular cases of linear parabolic equations 
with constant coefficients even the use of a three-point formula to 
approximate the first derivative, yielded unstable solutions (Reference 
47pp. 125 - 131). Though higher-order finite difference schemes can be 
employed, their stability and convergence may not be present, and due to 
the lack of availability of any comprehensive theoretical work for a 
non-linear system of equations with variable coefficients such as under 
discussion, their convergence and stability must be tested individually 
by actual numerical computation. 

For’the cases of ideal excitations, such as a large step of 
current and voltage considered here, an overall control of the truncation 
error is not possible, since higher-order derivations are either infinite 
or discontinuous. Step functions are represented in the’ discretized 
context as ramps since only non-zero step sizes can be used. The 
smallest mesh size can be determined by the usual limitations on the 
permissable range of the computer model. When the results for the cases 
of step functions were considered for different time-mesh sizes (within 
limits which guarantee stability), it has been observed that they are 
relatively insensitive to the variation of time step size (within limits) 
except for the points in the very vicinity of the discontinuity, where 
higher time derivitives exhibit large magnitudes. The considerable error 
introduced at these singular points may be restricted to very short time 


intervals by the use of very small time steps. A valuable feature 


consistently observed is that the error introduced in the neighbourhood 


of the singular points does not affect significantly the response at 


later times. To sum up, it may be said that the discretization error 


of the solution obtained by the formulation described in previous 
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sections, may be considered sufficiently small for any practical purposes 
even for step functions (numerically worst case) if the proper choice of 
time step-size is made. The choice of time step~size is governed by the 
abruptness of the variation of the external excitation, and LeSee treet 
on the time-variation of the unknowns, and is confirmed by testing on 
the actual numerical calculations for the evolution of the transients. 

6.8.2 Iteration Error 

The iterative procedure used at each instant of time, and 
described in the previous sections, are used to resolve the non- 
linearities of the problem by achieving consistency between the spatial 
distribution of the unknowns and the system of the discretized equations. 
When the estimated discrepancy between the computed distributions at the 
completion of an iterative cycle, at a given instant of time, and the 
‘true’ distributions, is less than a specified amount (determined by the 
conditions of convergence, and stability) the procedure restarts at the 
mext time instant.< This discrepancy represents the iteration error at a 
given time. The estimation of the stage, when overall consistency has 
been reached, can be made by either of the following criteria (depending 
upon the external sea neattion sy: 

(1) This estimation may be obtained by an appropriate 
extrapolation of the corrections for the distributions available from 
the previous cycles at the same instant of time. 

(2) A constant fraction of the maximum excursion (overall 
spatial points) of each unknown, at each instant of time, may be 
conveniently specified as the iteration error. 


(3) The estimation may be made by considering only the 


consistency of the terminal response at each instant of time which 
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which reflects sufficiently well the behaviour of internal distributions 
(except at singular points). 

(4) A more accurate evaluation of the consistency at a given 
time instant may be performed with the aid of a set of test relations 
derived from fundamental Eqs. (4.18) and (4.22). (The principle is the 
Same as was used in the steady-state case). Due to heavy computation 
load, this technique is not suitable to test consistency at each 
iteration cycle, but may be used to refine the accuracy of the solution 
Biya particular instant of time, for.which spatial, distribution of 
dependent quantities (such as particle current and net charge densities) 
are also required. 

6.8.3 Errors in the Initial Solution 

When the time dependent case is initiated from a steady-state 
condition, initial spatial distributions of unknowns must be furnished. 
These results are obtained from the solution for the steady-state case 
as described in Chapter V. The inaccuracy of these steady-state 
distributions represents a contribution to the numerical error of the 
transient solution. As such, it may generate unstable solutions, and 
therefore must be contained to sufficiently small values. It is the 
reason for taking all the precautionary steps we took to improve the 
accuracy of the steady-state solution. 

There is still another aspect to be considered. The analytical 
and discretized formulations adopted for both the steady-state case and 
the time-dependent cases are not 'compatible' to each other, in the 
Sense that the steady-state distributions feature truncation errors 
different from those introduced by the algorithm that generates transient 


Solutions. Ideally, if a transient solution is initiated from a steady- 
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state condition in the absence of an external excitation, the internal 
distributions and terminal properties should remain unchanged. Instead, 
as a consequence of the above mentioned incompatibility, a "drift" of 
the original steady-state distributions is observed throughout the 
fictitious transient solution. The drift continues until compatibility 
between the spatial distributions and the new discretization scheme is 
achieved. Although the magnitude of the total shift is relatively small 
fof the order of’ the truncation error), it is enough to yield incorrect 
results for the initial period of the transient for the dependant 
quantities (such as particle and displacement currents spatial 
distribution). For the later part of the transient, this incompatibility 
did not affect the results significantly. 
| As a general rule, we can say that when initial conditions 
incompatible with the time-dependent formulation are used, a preliminary 
computation phase in the absence of external conditions is required, 
before initiating the actual transient solution to achieve overall 
compatibility. In our own case we followed this rule to avoid the 
errors arising from 'incompatibility' of the initial solution. 
6.8.4 Other Sources of Numerical Error 

(1) Round-off error is always present in a numerical solution 
dealing with quantities represented by a finite number of digits. It is 
minimized by using double-precision arithmetic in numerical computations 
which carry sixteen decimal places as in the IBM 360/67 model on which 
we performed our computations. 


(2) Small differences between nearly equal numbers also 


contribute to the numerical error and are present in the expression for 


net charge density. 
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It may be stated that these two sources of error are not 
responsible for significant inaccuracies in the final results, since 
other contributions to the numerical errors (truncation and iteration 
errors) are dominant, 

6.8.5 Time Step, Iteration and Truncation Errors and 

Numerical Error Growth in the Time Domain 

A dominant factor, governing the rate of convergence at a given 
instant of time, is the magnitude of the time step at that instant. 
Larger steps are always responsible for slower convergence rates. The 
time step size, directly related to truncation error in the time domain, 
becomes a determining parameter very much affecting also the magnitude 
of the iteration error if the computations have to be completed within 
a certain amount of time. In such cases, iteration error may become 
dominant over the truncation error in the time domain, so that it is 
quite advantageous to select the time step size on the basis of the 
requirements dictated by a minimization of the iteration GLLor, ine 
total number of spatial points used also affects the convergence rate 
of the iterative procedure. This phenomenon is related to the 
propagation of the inconsistencies of the approximate distributions 
throughout the length of the device during the execution of an iterative 
cycle at a given time. A smaller number of spatial points with certain 
limits (chiefly dictated by the truncation error in the spatial domain) 
usually neers the convergence. In light of the above arguments a 
larger truncation error in the spatial domain is tolerated to decrease 
the number of spatial points, increase the convergence rate and 
consequently reduce the iteration error and thus save computation time. 


The abruptness of the doping profile is indirectly responsible for a 
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feouser iteration error im .the vicinity of the metallurgical interface. 
This is a consequence of the larger truncation error (in the spatial 
domain), occuring near the interface when we choose large spatial steps 
to increase the convergence. In actual practice, when this principle 
was applied to the structures considered, the results obtained were 
found to be satisfactory from any practical consideration. 

The accumulation, throughout the evolution of the transient 
solution, of the iteration error introduced at each instant of time is 
of great importance. When a constant fraction of the excursion of the 
unknowns within a time-step is specified as the iteration error, a 
transient solution featuring a monotonic pattern is affected by the same 
relative error if the single iteration errors are accumulative. In the 
case of diode No. 1, a considerably more favourable situation has been 
observed, The iteration error at a given time instant is partially 
compensated for by a larger excursion of the unknowns at the next instant 
of time, so that the resulting error is less than the sum of the 
iteration errors introduced at the previous instant. DeMari also 
observed this effect in his calculations on highly asymmetrical abrupt 
germanium p-n junctions for his simpler model and he attributes it to 
re encequence of the deviation from the "exact" condition of dynamic 
equilibrium at a given time caused by the iteration error, and he 
regarded it as an automatic feature inherent in the formulation adopted 
for the LSeera However, in the case of diode No. 2, iterative errors 
exhibited an accumulative behaviour, and thus endangered the stability 
at later parts of the transients, unless certain precautions yd iee a 
Moreover, due to the nature of the time variation of the independent 


variables, truncation errors in the time domain were found to be much 
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worse compared to the case of diode No. 1. Asa consequence, we were 
forced to choose a much smaller time step size. (of thé order.of 10) in 
the case of diode No. 2 compared to the case of diode No. 1. This 
allowed us to keep the truncation error in the time domain and the 


iteration error at a given instant of time well within limits of 


numerical instability. 


6.9 Summary 


An iterative method of the solution of the fundamental equations 
(Eqs. (4.18) - (4.23)) taken to govern the behaviour of the varactors 
under a time-dependent excitation has been presented in this chapter. 
The fundamental set of equations have been rearranged to a more 
convenient equivalent set of three non-linear partial differential 
equations (Eqs. 6.1, 6.4 and 6.5). Equations (6.4) and (6.5) are of the 
parabolic type, non-linear partial differential equations. Extrapolation 
from the simpler cases, treated by the theory of the numerical analysis 
lead to the choice of a generalised pure implicit discretization scheme; 
since restrictions on the discretization mesh inherent in explicit 
discretization scheme and sensitivity of generalised Crank-Nicholson 
discretization scheme to the numerical errors are not desirable. Gauss 
elimination and backward substitution method has been used to solve the 
system of linear algebraic equations generated as a result of the 
discretization. Details of the iterative procedure employed to resolve 
the various non-linearities inherent in the mathematical equations under 
consideration, has been presented. Then the numerical algorithm has 


—dri id the case when an 
been extended to cover the voltage-driven case al . 


equation relating the terminal voltage and current is known, Various 
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sources of numerical error have been outlined and briefly discussed. It 
has been postulated that truncation error in the time domain and 
@eeration error are mainly responsible for the overall accuracy of the 
results in the time domain. Therefore, the iteration and truncation 
errors have to be kept small. In view of the above requirement, a 
larger truncation error has been tolerated in ane spatial domain by 
choosing a smaller number of points into which the diode is divided. 
This enabled us to keep the time mesh size small without consuming too 
much computation time. It is pointed out that due to the nature of the 
time derivatives of the various quantities inside the device, diode No.2 
requires a smaller time-mesh size compared to diode No. 1 and hence 


diode No. 2 consumes more computation time. 
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CHAPTER VII 
RESULTS FOR STEADY-STATE CASE 


In this chapter, the computed results for both the internal and 
external properties of diode No. 1 and diode No. 2, when the forward 
steady-state terminal voltage of the diodes is specified, are presented 


and briefly discussed. 


fee introduction 

The mathematical equations which are taken to govern the 
behaviour of the diodes under steady~state conditions are outlined in 
Chapter IV. The numerical solution of these equations is discussed in 
detail in Chapter V. The internal distribution of the Tdesenient 
variables of the mathematical formulation (electrostatic potential, 
electron quasi-Fermi level and hole quasi-Fermi level) and of the 
dependent variables (electron density, hole density, electron current 
density and hole current density) are useful in explaining the terminal 
behaviour. They are thus indirectly helpful in device design. The 
mobility and the generation-recombination rate distribution inside the 
device expose the validity of the physical model in which they are 
assumed constant. Moreover, the internal distributions of the hole 
density, the electron density and electrostatic potential are required 
for time-dependent solutions initial values, as already discussed in 
Chapter VI. The static characteristics of varactors and of step-recovery 
diodes, such as terminal current and incremental capacitance as a 


function of terminal voltage are very significant, though not sufficient 
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be chayacterize them.! All the results are displayed in normalised 


units, unless otherwise specified (see Table 4.1). 


7.2 Results For Diode No. 1 

Detailed results for the internal distributions for both the 
independent and dependent variables are displayed in this section for 
Beiminadevoltages of 0,0, 0.7, 0.9 and 1.0 volts. These results were 
stored on magnetic tape and were plotted directly by means of the 
Calcomp plotter attached to the IBM 360/67 system. The results displayed 
age electron density distribution, hole density distribution, net doping 
profile, electrostatic potential (E.S. Potential) distribution, electron 
quasi-Fermi level (E.Q.F.L.) distribution, hole quasi-Fermi level 
C70. F.L.) distribution, electron current distribution, hole current 
distribution and generation-recombination rate distribution throughout 
the interior of the device. 

Pieures (7 ala), 7.1b), (7.le) and (/eld) eive the electron) and 
hole density distributions along with the net doping concentration 
distribution. These results reaffirm the inequality between the hole- 
electron density product and externally applied voltage es 
pes) n(x) < eva, This is in agreement with aunties: and DeMari's 
recent numerical results for highly asymmetrical p-n junctions. The 
concentrations of minority-carriers near the external contacts is of 
interest as we have assumed ohmic contacts in our model. The 
concentration of minority-carriers increases as forward-bias is 
increased. Even at 1.0 volt forward-bias, the accumulation of minority- 
carriers near the ohmic contact at the substrate is relatively 


insignificant. On the other hand, at the ohmic contact at the end of 
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Fig. 7.la. Spatial distribution of net doping and mobile carrier densities 


for diode No. 1 under steady-state conditions for Va = 0.0 volt. 
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DIODE No. 1 
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SOLE 


rn iy 12 | 18 | 24 30 
DISTANCE (microns) 


lb. Spatial distribution of net doping and mobile carrier densities 


for diode No. 1 under steady-state conditions for Va = 0.7 volt. 
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Fig. 7.lc. Spatial distribution of net doping and mobile carrier densities 


for diode No, 1 under steady-state conditions for Va = 0.9 volt. 
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Fig. 7.ld. Spatial distribution of net doping and mobile carrier densities 


for diode No, 1 under steady-state conditions for es 1 Onvols. 
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the narrow p-region, there is a greater accumulation of the minority- 
eecriers (710 times more than on then side). But thie concentration 
of minority-carriers is small enough compared to the concentration of 
majority-carriers to significantly alter the charge neutrality near the 
external contacts on both the p- and n-side. It may be pointed out that 
at very high injection levels, when the minority-carrier concentrations 
become ae oe to the majority-carrier concentrations, neutrality will 
be absent in a narrow region near the external contacts. Since the 
assumption of ohmic contacts will prevent the conductivity of such a 
region from being significantly modulated, a large current will cause a 
high voltage drop or, in equivalent terms, a strong electric field near 
Premcontacts: to’ be. built up. This electric field will, in.turn; tend 
to expel the majority-carriers and attract the minority-carriers, thus 
resulting in a local alteration of neutrality balance. Conductivity 
modulation effects due to the minority-carrier injection, are absent till 
0.9 V and start showing their presence for 1.0 Volt forward-bias. The 
narrow region near the metallurgical junction shows accumulation of the 
minority-carrier densities greater than the majority-carriers densities. 

Fieuvese(1.2a)e) (7.26), (/42e) and (7.2d) display the Es. 
potential, pesO.F.L. and %0.F.L. distribution throughout) the interior of 
the device. These illustrations show that the assumption of constant 
Fermi levels in the transition region is a valid assumption for moderate 
forward-bias levels. Due to the injection of carriers across the 
junction, the abruptness of the E.S. Potential distribution at the 
metallurgical junction decreases gradually as the forward bias is 
increased. The E.S. Potential distribution for 1.0 volt forward bias 


(Fig. 7.2d) shows a small slope in the substrate. It shows that the 
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Fig, 7.2a. Spatial distribution of electrostatic potential and hole and 
electron quasi-Fermi levels for diode No. 1 under steady-state 


conditions for Vag = 0/0 Wolt. 
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Fig, 7.2d. Spatial distribution of electrostatic potential and hole and 
electron quasi-Fermi levels for diode No. 1 under steady-state 
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presence of the minority~carriers in the substrate has become 
significant enough to cause a slight shift in the neutrality condition 
present before. 

Figures (7.3a), (7.3b) and (7.3c) show the net generation- 
recombination rate for 0.7, 0.9 and 1.0 volts forward bias, respectively. 
Initially, most of the recombination is confined near the metallurgical 
junetion (Fig. 7.3a), but as the forward bias is increased, recombination 
in the quase-neutral region of the diode becomes more and more important. 
For 1.0 volts forward bias, (Fig. 7.3c) shows that recombination is 
significant not only in the transition layer, but also in the quasi- 
neutral regions of the diodes, 

Figures (7.4a), (7.4b), (7.4c) and (7.4d) show the mobility 
fase DUution thromghout} the interior of the device for 0.0, 0.7, 0.9 and 
1.0 volts forward bias. These figures show that the assumption of 
constant mobilities (as assumed in most of the theoretical work) is 
hichly inaccuvate in the vicinity of the metallurgfeal junction, where 
mars a strome function) of position due to the presence of the elettric 
field and to the variation of net doping density. The electron mobility 
is about 1.5 times greater than the hole mobility, as expected. 

Figures (775a)}, (7.5b) and (7.5c) show the spatial distribution 
of electron and hole parecne” for 0.7, 0.9, and 1.0 volts forward—bias, 
This shows an interesting feature in that, though the doping density of 
the p-side is higher (5 X ihe than that of the n-side (4 X 107°) , the 


electron current is the dominant current flowing in the substrate and 


ee ee ee 


: For simplicity, in the following text, the term "current" will 


be consistently used to denote current density. 
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Fig. 7.3a. Spatial distribution of net recombination rate for diode 


No. 1 under steady-state conditions for V, = 0.7 volt. 
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Fig. 7,.3b. Spatial distribution of net recombination rate for diode 


No. 1 under steady-state conditions for Vy = 0.9 volt. 
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Fig. 7.3c. Spatial distribution of net recombination rate for diode 


No. l under steady—state conditions for Vy = 1.0 volt. 
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fio, j4ay. Spatial distribution of carrier mobilities for diode No. l 


under steady-state conditions for Va = 0.0 volt. 
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ime. ¥.4b. Spatial distribution of carrier mobilities for diode No. l 


under steady-state conditions for V, = 0./ volt. 
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Fig. 7.4c. Spatial distribution of carrier mobilities for diode No. l 


under steady-state conditions for Vy = 0.9 volt. 
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Fig. 7.4d. Spatial distribution of carrier mobilities for diode No. l 


under steady~state conditions for Vy = Pe Oovolts, 
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Fig. 7.5a. Spatial distribution of electron and hole currents for 


diode No. 1 under steady-state conditions for V, = 0.7 volt. 
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Fig. 7.5b. Spatial distribution of electron and hole currents for 


diode No. 1 under steady-state conditions for V, = 0.9 volt. 


r, en agoic 


i, : LD ‘i - ; Hi * NWor, we iRicad a af home nace 
, 


j : vi ; 7 nf : o n — | ari 


vee age: 5. Sa aaah a i 


22 


10% ne 
{ Cileedent samen 
DIODE Nou! 
Vy 21.0 VOLT 
as ELECTRON-/ 
WS) 
— 
tis 
Ly MOLE 
OZ 
CE 
=) 
U 
107 
Z 
O = 
= 
YU 
Zz 
=) 
10° seas (eres eerie eae Poe ee) eh LY 


mee 4 12 18 24 
DISTANCE ( microns ) 


Fig. 7.5c. Spatial distribution of electron and hole currents for 


diode No. 1 under steady-state conditions for Vag = 1.0 volt. 
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makes a significant contribution to the total current in the p-side. 
The hole current in the p-region makes a decreasing contribution to the 
total current as the forward-bias is increased; it contributes ol 772 of 
the total current in the p-region for 0.7 V bias compared to 72.5% and 
woe. tor 0.9°V and 1.0 V forward bias, coe CHEN The contribution 
of the electron current to the total current in the p-region is 
increasing as the forward bias is increased. This may be explained in 
the following way: In our assumed model, the ratio of hole current to 


ebectron current can be written as 


wi) Jp diftusion (x)= J dritt. ) 
eRe a eis at Pha eer ns ee has 


Pehl sk) 
Jeo me rdcusion (x) oJ dere tks) 
n n n 
J_ (x) 
In the case under discussion, 7 Ty decreases in the p-region. As the 
. 


forward bias is increased, the electric field in the p-region also 
imereases, which gives rise to a large increase in the hole drift current 
Moles are majority carriers, and therefore their concentration is very 
high and a small change in the electric field will cause a large change 
in the drift current). It seems that as the forward-bias is increased, 
the increment change a the hole drift current is greater, compared to 
that of the hole diffusion current. ay a consequence of the above 
reasoning, we see that the hole current contribution to the total current 
will decrease in the p-region as the forward-bias is increased, The 
argument that the electron mobility is greater than the hole mobility is 
not scelevant to this.case, as their ratio is almost constant in the p- 
region for the biases considered (except in the transition region). The 


contribution of the hole current to the total current in the quasi- 
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meutral n-region increases slightly as the forward bias is increased 
(0.8, 1.27, 1.81% near the external contact at 0.7, 0.9 and 1.0 Volts, 
respectively). In this region, the diffusion component is increased 
more than the drift component as the terminal voltage is increased, 
piace the electric field is either absent or has a very small value. 

The terminal properties of the diode, namely the terminal 
current and the incremental capacitance as a function of the diode 
terminal voltage, are plotted in Figs. (7.6) and (7.7). The incremental 
capacitance starts showing a large variation at a forward-bias of 0.7 
volts, and increases very rapidly to a value approximately 10° times the 
value at zero volts. The value of the series resistance at zero volts 
bias calculated from the definition given in Eq. (5.32) is 0.81 ohms. 


: : 1 
Cut-off frequency, the static figure of merit, is defined as follows 


fo =—~—— BAC ioe) 

oa Jeg a Oe 
sv jv 
where 

f = cut-off frequency at the specified voltage Cv)" 

reecy. 

R = diode series resistance at the specified voltage (v). 
SV 

C. = diode junction capacitance at the specified 
ap 


voltage (v). 


The cut-off frequency for diode No. 1 at zero volts bias, calculated 
mom iq. (/42), 16.261. GHz. Ltyis to be noted that the above calculation 


of the cut-off frequency did not include the package capacitance and the 


eoptact resistance. 
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Fig. 7.6, Terminal current of diode No. 1 as a function of terminal 


voltage (steady-state conditions). 
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Fig. 7.7. Incremental capacitance of diode No. 1 as a function of 


terminal voltage (steady-state conditions). 
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feo Results For Diode No! “2 


The results displayed are hole and electron density 
distributions, electrostatic potential, hole and electron quasi-Fermi 
level distributions, hole and electron current distributions, hole and 
electron mobility distributions and generation-recombination rate 
distribution. These quantities are displayed for 0.0, 0.5, 0.8 and 1.0 
weres. Figures (7.84), (7.8b), (7.8c) and (7.8d) give the hole and 
electron density distribution along with the net doping concentration 
distribution throughout the interior of the device. These results 
affirm the inequality between the electron-hole density product and the 
externally applied voltage ee Ds) nay) < eva, which was earlier 
Pointed out by Gumal”> for simple p-n junctions, and was reaffirmed by 
our results on diode No. 1 as mentioned in the previous section. This 
result shows that in our at = i= n’ structure, the electron and hole 
density distributions obey the same law as for simple p-n junctions. 

The hole and electron concentrations show that a great number of the 
injected minority carriers are stored in the middle n-region. 
Consequently the concentration of minority carriers in the quasi~neutral 
n’ region is far less (by a factor of more then ten) than the 
concentration in the middle high resistivity layer. This may be 
explained with the help of Fig. (7.9), which displays the electric field 
in the neighbourhood of the middle n-region, when the diode is in thermal 
equilibrium. On both edges of the n-region, there is an appreciable 


a P A =e ’ 
negative electric field present. At the p - nm junction the field is 


ge ; 
about 5 times higher then that at the n-n edge. The electric field 


at the center of the n-region is zero. Due to the presence of this 


electric field, the flow of injected carriers into the n and p_ regions 
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Fig. 7.8a. Spatial distribution of net doping and mobile carrier 
densities for diode No. 2 under steady-state conditions for 


Va = 0.0 vole, 
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Fig. 7.8b. Spatial distribution of net doping and mobile carrier 


densities for diode No. 2 under steady-state conditions for 


Va = 0.5 volt.. 
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Fig. 7.8c, Spatial distribution of net doping and mobile carrier 


densities for diode No. 2 under steady-state conditions for 


Va = 0.8 volt, 
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Pig. 7.8d. Spatial distribution of net doping and mobile carrier 


densities for diode No. 2 under steady-state conditions for 
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is retarded, resulting in an accumulation of injected charge in and 
ecound themeregion. This is a desirable effect since it tends to 
localize the injected charge in the vicinity of the {unceLlon, VGhas 
allowing it to be recovered in a shorter time when the diode is driven 
into reverse-bias (discussed in detail in Chapter IX). It may be 
pointed out that due to the presence of the electric field in a quite 
wide region, more energy will be dissipated in the diode during the 
transit times and will contribute another loss component to the total 
loss. The accumulation of the charge in the center n-region causes 
conductivity modulation, since the injected carrier concentration is 
much higher (of the order of 10° times for 1.0 volts) than the doping 
density (1.0 X Loy J ana 

The concentration of minority carriers near the external 
contacts exhibits a behaviour similiar to that observed in diode No. l. 
The carrier concentration near the contacts increases as the terminal 
voltage is increased and is much greater (of the order of 107 times) on 
the p -side than on the ede The concentration is observed to be 
less—than that in diode No. 1, for the same current. It may be 
attributed to the presence of a large negative electric field, spread 
over a larger distance. 

Figures (7.10a), (7.10b), (7.10c) and (7.10d) display the 
distribution of E.S. Potential, E.Q.F.L. and H.Q.F.L. throughout the 
interior of the device. The figure shows that the E.S. Potential 
distribution has a large slope (hence electric field) aE the junctions 
between the regions of different doping. This slope (hence electric 
field) decreases as the forward bias is increased. Another important 


feature observed is the zero slope of the E.S. Potential (zero electric 
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Fig. 7.10a. Spatial distribution of electrostatic potential and electron 


and hole quasi~Fermi levels for diode No. 2 under steady-state 


conditions for Va = 0,0 volt. 
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Fig. 7.10b. Spatial distribution of electrostatic potential and electron 


and hole quasi-Fermi levels for diode No, 2 under steady-state 


conditions for Va > On 5S vob. 
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Pie! 7 S106 Spatial distribution of electrostatic potential and electron 


and hole quasi-Fermi levels for diode No. 2 under steady-state 


conditions for Va = 0.8 volt. 
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Fig. 7.10d. Spatial distribution of electrostatic potential and electron 


and hole quasi-Fermi levels for diode No, 2 under steady-state 


conditions for ee PO vor‘, 
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field) in the vicinity of the center of the n-region. This charge- 
neutral region expands towards the enon as the forward bias is 
increased. If we look at Figs. (7.8a), (7.8b), (7.8c) and (7.8d), we 
see that the mobile charge density of both the hole and electrons 
injected into the n-region is equal (the doping concentration of the 
n-region is very low compared to that of the injected carriers), hence 
charge-neutrality prevails. The expansion of the central charge-neutral 
region towards the peer ceion shows that the injection of the electrons 
becomes more and more dominant with increasing forward—bias. This 
argument is supported by the results for the electron and hole current 
distributions. At 1.0 volts forward-bias (Fig. 7.10d), a small positive 
electric field is present throughout the n' +region. Hole quasi-Fermi 
level also shows an appreciable slope in this region, which indicates 
the heavy injection of holes. 

Thece resules ofthe, B.S. Potential, the £0. F.U. end t.Oa. &. 
distributions expose various inconsistencies in the analytical nodet-: ec 
The analytical model divides the diode into charge-neutral and space- 
charge regions, where space-charge is confined to the narrow layers at 
the junctions between regions of different doping. The E.Q.F.L. and 
H.Q.F.L. are assumed constant throughout both the space-charge regions. 
This model seems to fail at high forward-bias, where the electric field 
in both the . 2 and p'-regions is no longer a negligible quantity, and 
the E.Q.F.L..and H.Q.F.L. are no longer constant in the vicinity of the 
junctions between regions of different doping. 


The distribution of the generation-recombination rate throughout 


the interior of the device is plotted in Figs. (7.lla), Cia 1 iboLeand 


(7.1lc). For low forward bias most of the recombination takes place in 
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mie, 7.1ia, Spatial distribution of net recombination rate for 


diode No. 2 under steady-state conditions for Va = 0.5 volt. 
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Fig. 7.1llib. Spatial distribution of net recombination rate for 


diode No. 2 under steady-state conditions for Vy 


= 0.3 volt. 
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Fig. /,lle. Spatial distribution of net recombination rate for 


diode No. 2 under steady-state conditions for Va = 1.0 volt. 
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the n-region. As the forward bias is increased, the Deeeon is Stal 
the region having greater recombination, but the region of greater 
recombination expands further into the ee yack Figures (/.12a), 
eb), (7.120) and .(7.12d) display the hole and electron mobility 
distributions. The mobilities of both the hole and electron are about 
nine times greater in the n-region than in the quasi-neutral oo and ae 
regions. This is due to the fact that the n~region is very lightly 
doped compared to the ye and Sone, Moreover, at low forward- 
biases, the high electric field at the a junction causes the 
mobilities to decrease appreciably in the vicinity of the junction. 

Figures) (7,13a), (7.13b) and (7.13c) show the hole and electron 
eucrent distribution for 0.5, 0.8 and 1.0 volts respectively. Initially, 
the hole and electron currents dominate the total current in the ou and 
n’ regions, respectively. As the forward bias is increased, the 
electron current exhibits a marked increase in all three regions, until 
it becomes the major contributor of the total current throughout the 
@uode (more than 90Z at 1.0 volt). The contribution of the hole current 
in the py venton for 0.5 V forward-bias is about 96%, decreasing to 
about 76% at 0.8 V. At high current densities (1.0 V forward-bias), it 
reduces to 5.16%. This marked decrease of the hole current contribution 
to the total current may be explained in the same way as explained in 
Section (7.2) for diode No. 1. The diffusion and drift components of 
the hole current are opposite in direction. As the forward bias is 
increased, the diffusion component increases, but as it appears that the 
drift components have a larger increment (holes are majority carriers in 
the pe and a small electric field is enough to give a large electric 


current) than the diffusion component, particularly at high forward-bias, 
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mie, 7.12a. Spatial distribution of carrier mobilities for diode No. 2 


under steady-state conditions for Va = 0.0 volt. 
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Fig. 7.12b. Spatial distribution of carrier mobilities for diode No, 2 


under steady-state conditions for Vy = 0.5 volt. 
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Fig. 7.12c. Spatial distribution of carrier mobilities for diode No, 2 


under steady-state conditions for Va = 0.8 volt. 
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Fig, 7.12d. Spatial distribution of carrier mobilities for diode No. 2 


under steady-state conditions for Va = i VOLS 
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Fig. 7.13a. Spatial distribution of electron and hole currents for 


diode No. 2 under steady-state conditions for Va = 0.5 volt. 
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ig. 7.13b. Spatial distribution of electron and hole currents for 


diode No. 2 under steady-state conditions for V, = 0.8 volt. 
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Hig. 7.13c. Spatial distribution of electron and hole currents for 


diode No. 2 under steady-state conditions for Va = 1.0 volt. 
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when the electric field in the pl wresion becomes more and more signifi- 
cant. It results in less and less hole current contribution to the 
total current as the forward-bias is increased. At low and moderate 
forward-biases, most of the decrease, in the hole and electron currents 
takes place in the central n-region. It may be attributed to the 
presence of high negative electric fields present at both ends of the 
n-region. At high forward-bias, when the space-charge regions at the 
junctions of the regions of different doping are partly neutralized, no 
such abrupt variation takes place at the ends of the n-layers. 

The eee properties of the diode, namely the terminal 
current and the incremental capacitance as a function of the diode 
terminal voltage, are plotted in Figs. (7.14) and (7.15). ~ The value of 
the aa frequency at zero bias calculated fromiiq. Cr .2):, is 
52.6 GHz. It is to be noted that this value of the cut-off frequency is 
the upper limit that may be expected, since we did not include ‘the 


package capacitance and contact resistance. 


7.4 Synopsis 


The internal and external behaviour of diode No. 1 and diode 
No. 2, for steady-state voltage excitation, has been presented. The 
results, displayed in graphical form as a function of position and of 
applied terminal voltage, give a thorough understanding of the devices' 


behaviour. A brief discussion of the results has indicated a number of 


interesting features and has demonstrated the usefulness of this 


computer simulation method in making a wealth of information available. 
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Fig. 7.14, Terminal current of diode No. 2.49) a Lunctionsof terminal 


voltage (steady-state conditions). 
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CHAPTER VIIL 
LARGE-SIGNAL "TURN-ON" RESPONSE 


In this chapter, the computed results for both the internal 
and external properties of diode No. 1 and diode No. 2, when "turned-on" 


by a large step of current, are presented and briefly discussed. 


8.1 Introduction 

When the bias of a diode is switched from either the reverse or 
the equilibrium state to a forward state, the behaviour of the diode 
under the condition is termed the "turn-on" response. From the point of 
view of charge-storage behaviour, this phenomenon is important. Unless 
a sufficient amount of charge is injected and stored during the forward 
Pecerot the a.c,-cycle, efficient use of the device as a step~recovery 
diode, inspite of a very fast recovery process, is impaired. 

It is the purpose of this chapter to examine in detail and to 
explain quantitatively the forward transient behaviour of diode No. l 
and diode No. 2, which are taken to represent the two distinct kinds of 


varactor diodes as discussed in Chapter II. 


A number of workers have studied this problem for a highly 


asymmetrical step junction diode. The reported characteristics serve as 


° . . ° WwW 
the guideline in this chapter, though our case is "doubly-injecting 
Priecoth carriers are to be considered . The two important resulte whien 
are reported in the literature may be summarised as follows: 


(1) The initial condition of the diode, either in the reverse 


state or in thermal equilibrium, has little effect on the observed 
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forward transient behaviour. But, if the diode is initially biased 
slightly in the forward direction, the maximum deviation between the 
initial value and the final value (voltage or current) can be greatly 
reduced; 

(2) The diode forward transient is considered to consist of 
two components; one behaves in a manner similar to an inductive circuit, 
and the second in a manner similar to a capacitive circuit. For large 
currents, the inductive character dominates, while for small currents, it 
is the capacitive behaviour that dominates. Between these two extremes 
the transient exhibits an oscillatory behaviour. 

The case discussed in this chapter is the behaviour of the 
diodes when a large step of current is suddenly applied across the diode 
terminals. All the results are plotted in normalised units, unless 


otherwise specified. (See Table 4.1) 


6.2 Results for Diode No. 1 


A current step having a value 


0.16 X 10° E70 


ea 
il 


is applied directly to the terminals of diode No. l, initially in thermal 
equilibrium condition. This value of current corresponds to 1.0 volts 
forward bias under steady-state conditions. The overall terminal voltage 
response is shown in Fig. (8.1). Figure (8.2) displays the response for 
the initial part of the transient. It is governed by the ohmic drop 


build-up and the charging of the depletion layer capacitance. The build- 
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up of the ohmic drop takes place in a time of the order of two times the 
dielectric relaxation time (1.1 X Fanos sec. for the substrate having 
0.4 X ion ohm—cm, resistivity). Since the resistivity is so low, the 
ohmic-drop build-up is very small and very fast. Therefore, its 
influence on the overall transient is insignificant. 

The overall transient may be taken to consist of two components: 
the junction voltage transient and the diode resistance transient. The 
former increases with time and the latter decreases with time. In our 
case, charging of the depletion layer capacitance in the initial stages 
and of the diffusion capacitance at later times dominate the overall 
transient response. The conductivity modulation which results in a 
decrease in the diode resistance has little effect even at latter times 
when a large number of the minority-carriers have been injected across 
Eqesjunction, It is due to the facts. that tthis diode is very heavily 
doped and for any conductivity modulation to take place, the injected 
carrier density is to be comparable to the majority-carrier 


concentrations. 


Figure (8.3) displays the electric field distribution for 
various time instants during the evolution of the transient. Figures 
(8.4) and (8.5) display the electron and hole density distributions for 
several time instants (the same as chosen for the display of the electric 
field) during the evolution of the transient. As is obvious from Fig. 
(8.3), the displacement current which is initially the major current | 
aoe throughout the diode, becomes negligible in the quasi~-neutral 
regions within a short time of the order of 2 X Tent Ct is the 
dielectric relaxation time of the A recton)s but is still significant 


in the transition layer. However, it gradually decreases even in the 
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Bic. 2.3 Evolution of the electric field distribution in diode No. 1 during 


the turn-on transient (values of time instants indicated in Fig. 8.1). 
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Eransition layer towards its steady~state value of zero as the transient 
evolves in time. At this stage the transient is governed by the motion 
of the minority-carriers, and hence becomes slow in time. The diode 
will approach the steady-state conditions with a time constant of the 
order of the transit time of the holes in the substrate (=1.25 X Lopeee):. 

Figure (8.6) displays the ratio of the hole charge injected 
across the junction into the n-region to the holes stored in the n-region 
under the steady-state conditions for the same terminal current. It is 
to be noted that less than 5% of the charge is injected in a time of 
igs pa oa. This time is comparable to the time period of a 4 GHz. signal. 
Moreover, this ratio will approach its final value of 1.0 with a time 


constant of the order of the transit time of the holes in the n-region 


1.5: X 107° sec.) 


Peo nesults for Diode No. 2 


A current step having a value 


=s0,/9 4 10? te) 


(4 
| 


is applied directly to the diode. The diode is assumed to be in the 
thermal equilibrium before the current step is applied. The value of 
the current corresponds to 0.915 volts forward-bias under the steady- 


state conditions. The overall voltage response is displayed in Fig. 


(8.7). There is a overshoot present in the voltage response. The 


-11 ; 
initial part of the transient (t < 1.1 X 10 sec.) mainly represents 


the build-up of ohmic drop in the quasi-neutral regions of the diode and 
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the charging of the transition region capacitance, This phenomenon is 
ruled by dielectric relaxation properties of the high-resistivity region, 
and occurs within a period approximately equal to the dielectric 
relaxation time, Tw as a consequence of the dielectric relaxation 


times of the surrounding low resistivity regions ail and a are very 


Ee 


small compared to ae Ac-later Himes Cparticularly it 70.24 et 00 
the conductivity modulation effects, comparatively small in the initial 
part of the response, are dominant, so that the diode terminal voltage 
decays to its final steady-state value. This decay is mainly governed 
by the conductivity modulation of the high-resistivity n-region, since 
in the substrate, the effect is far less because of its initial lower- 
resistivity and the fact that the injected carriers have to pass through 
the n-region before they reach the n’-region. Thus the time constant of 
the decay is essentially given by the transit time of the holes in the 
n-region (6.0 X Kone ae Note that the transit time for electrons is 
smaller due to their higher mobility. (=2 times that of holes) 

Figures (8.9) and (8.10) display the hole and electron densities 
during the evolution of the transient. In Figs. (8.lla) and (8.1lb), the 
electric field as a function of distance is plotted for various instants 
of time. Figures (8.12a) and (8.12b) give the displacement current as a 
function of distance for various instants of time during the evolution 
The electron drift current component in the mn -region 


Of the transient. 


as a function of distance is displayed in Fig. (8.13) for various 


* This estimation is made from the expression T,. TT 


W is the estimated width of the n-region and Dp is the average 


diffusion constant of holes in this region. 
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Lay (3) 
instants of time during the initial part of the transient. As is evident 
from Eqs. (8.12a) and (8.12b), displacement current features a basically 
different pattern in the n-region when compared to diode No. 1. The 
initial decrease of the displacement current, which corresponds to the 
build-up of the ohmic potential drop throughout the diode, is followed 
by an inversion of sign in the n-region and a subseguent decrease in the 
magnitude towards its steady-state zero value. The reversal of the 
direction of the displacement current flow is a consequence of the 
decrease of the voltage drop in a region whose conductivity is being 
modulated. The modulation initiates in the n-region at the edge of the 
transition region of the oe junction and gradually spreads towards the 
Beorecion, The conductivity is well apparent in the spatial distributions 
of the mobile carrier densities (fies. 6.9 and 3,10) ) 9 Due wte smaller 
dielectric relaxation time, the displacement current, dominant in the 
initial stages, decreases markedly in the n’-region, though Tt temstivy 
dominant in the n-region. As shown in Fig. (8.13), the electron drift 

P ip , ves 
current replaces the displacement current in the n -region. This is the 
consequence of the small electric field built up in the n'-region due to 
the flow of the current. Similiar behaviour is exhibited by the ace 
region. The electric field in the n-region exhibits an interesting 


behaviour (Figs. 8.lla and 8.1lb). Initially the electric field is 


either negative or zero in this region. As the transient evolves in 


se 
time, the electric field reverses its sign at the n-n edge and becomes 


+- 
positive. This positive electric field region spreads towards the p - 


—5 ; 
region edge. At later time (t ~ 0.28 X 10 ~) the whole n-region has 


: + me ONE 
electric field and even the transition region (p -n) exhibits a large 


positive electric field. But as the transient proceeds further in time, 
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ate : . 
the electric field at the p -n junction edge starts becoming less 


positive. At this point in the time domain, the decay of the terminal 


voltage towards its final steady-state value starts. 


8.4 Conclusions 

Diode No. 1 and diode No. 2 exhibit basically different terminal 
behaviour. The overall transient of diode No. 1 behaves some what 
similar to an RC circuit. Initial large conductivity of the semi- 
conductor material on both sides of the p-n junction prevents 
conductivity modulation effects from becoming significant at the current 
density to which the diode is subjected. It is also shown that less 
than 54 of the holes stored for this current density under the steady- 
Seate conditions is injected in 2.2 X ee seconds after the transient 
is initiated. The displacement current is dominant in the initial part 


of the transient throughout the diode. In a time of the order of two 


16 


times the dielectric relaxation time of the substrate (+1.1 X 10 ‘sec.), 


the particle drift current takes over from the displacement current in 
the quasi-neutral bulk regions, but is still significant in the 
transition region. At later times, the minority-carrier injected across 
the junction become significant and the transient is governed by their 
flow. The steady-state value will thus be reached with a time constant 
of the order of the transit time of the holes in the n-region 
(led ae Se 

The overall transient of diode No. 2 behaves somewhat similiar 
| much greater than the steady- 


to an RL circuit, the initial voltage being 


state voltage. The conductivity modulation effects in the high- 


resistivity n-region are mainly responsible for this behaviour. The 


\ 


\. sein : . >. home x3 ost Xb ietive oid) 


betel: sieck 2 oghd. 2 arte ol th. 


~ ceubearoe Lonfth of) af-aien 


aoteuda 

—: - 

. oo 

ofl eholi bon I cof SbOke 


.~Stwosls. 08 
vie djed ne leatisiam., 


eforiis sobteiubom Caza 

2) shoth oft dabshr ot “2. 
ay 

ssa eeloni ott Woe 

¥ Aa 4 


tosis i anoti lage 


ATqnth oft .bedakae 


\ 


‘sotgee rds Joslengya ge 
ie 


cc +3qngisa abs tootelbe 4 
e 7) 

: f 

esus Jog Eo J2i7b slakiy 
‘ ies 
,anoiged) hud thsssda-taee oP 
ent 

. (+ sétel JA .oolges wo) 3 Tee 
; a 4 

biogiea smoasd hu 9 

> -: 


[hy 62 ak Vadim He mae: 


iis 


at 


Spoke: Fm oS pe jet Siexeve. at her 
rth i " : na a : 
) ete ; ne, © 


a we & 

ary - pare Ag ne 

deapttiov Deaths acne? 7 a 

f is TY | i) oe 

‘’ aot ie hes os ahs ; wala - 
vottehibon yjiedaahioag $e% > ogestov tap 
aaeeaet et. ee 


scvieladts ilannlasiias Sata ileal 


* on 7 ~ = 
f aa cy r EAS Tee Vi3 2 ee ‘ 
oii adi vena 
P 4 bd ans = her a” . eee 
4 - le al ; ey ot, sh 4 af} i : 
r i Ua , See La iY . f oT 


EE cil Vit 


3) 
large forward current to which the diode is subjected, causes a high 
positive electric field to build up in the n-region. It later spreads 

ct : : Dane . 
to even the p -—n junction transition layer. At later times 


ae x 1077 


sec.), due to significant increases in the conductivity at 
Cie e * ° e e ° 2 

the p -n junction edge, the positive electric field starts decreasing. 

At this point, the terminal voltage starts decaying towards its final 


value. This decay will be mainly governed by the flow of the mobile 


Carriers in the n-region. 


Oh 
= t oa ,(.28e oF Ae 7 
1 d a 
i ' ,sghe sobktoayes 2 


Avdide of so Tawe Ree “4 


Wie of fad atyioels 
uD 


ee 
2013200), i - G Sia ae 


ease m3 ext3 sHteg Rage ac ' 
7. 


as 


6 1i:2¥ VEISD eit oat 


/ weige2a- elt ae 


Ra 
' 
SS 
/ 
4! 
! * ae 


176 


CHAPTER IX 
LARGE-SIGNAL "TURN-OFF" RESPONSE 


In this chapter, the computed results for both the internal and 
external behaviour of diode No. 1 and diode No. 2, when the device is 
"turned-off" with a current-limiting series resistance, are presented 


and briefly discussed. 


we4 Introduction 
When a diode has been carrying current in the forward direction 
and a reverse-biasing potential is suddenly applied through a series 
resistance, the behaviour of the diode is called the "turn — off" 
response. We are interested in this phenomenon to study the recovery of 
the charge injected during the initial forward-biased condition. The 
purpose of the series resistance is to limit the initial reverse current 
which will start flowing instantly after the reverse-biasing potential 
is applied. The time required by the diode to recover its stored charge 
determines to a great degree the upper limit of the frequency at which 
a diode can be operated in the step-recovery mode. 
A number of workers have studied this problem for a highly 
; 9-13, 16 


asymmetrical step p-n junction, Their analysis has been mainly 


based upon the solution of a greatly simplified diffusion equation, 
which does not take into consideration the presence of the electric 
field. Muto et Bie expanded the analysis for a graded-base diode and 
found it to have a faster response than a comparable uniform-base diode. 


feoredy presented a second-order analysis taking into account various 
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recovery mechanisms which were neglected in the previous works mentioned 
above. In the analysis of Benda et al.’’; a theory is developed for a 
p-S-—n rectifier. The center high resistivity layer in diode No. 2 
isivery thin. Therefore, the theory of Benda et al!’ is not applicable 
as such and has to be modified. 

In this chapter, the results for a "turn-off" transient for 
both diode No. 1 and diode No. 2 are presented. These results are based 
on our physical and mathematical model as discussed in Chapter III and 
Chapter IV. The details of the numerical solution are discussed in 


Chapter VI. All the quantities are displayed in normalised units, unless 


otherwise specified. (See Table 4.1) 


9.2 Results for Diode No. l 
The parameters chosen for the switching circuit of Fig. (9.1) 


for diode No. 1 are as follows: 


Ve = —20 volts 
Ve = +1 volt 
Rel OG 


Figure (9.3) displays the diode terminal voltage for short times 


(<10 !! sec.) after the diode is subjected to the "turn-off" transient. 


Ve + Ve CE=0) 


us RPG FN, ON F 
At the time instant t = 0, a current I(t=0) = 2 is 
forced through the device, since its terminal voltage remains constant 


from t=O to t=0'” As soon as the current starts flowing, the build-up 


of ohmic drop in the quasi-neutral p and n regions becomes the dominant 


effect. The build-up of ohmic drop effectively takes place in a time 


constant of the order of the dielectric relaxation time of the n—-region 
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Fig. 9.1. Circuit model used for the study of the large-signal 


turn-off response of the diodes. 


Fig. 9.2. Switching waveform applied to the diodes. 
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-16 
ct, X 10 sec. for an average doping density of ee Cie The 


approximate value of the ohmic drop in the n-region may be estimated as 


follows: 
_ P& iS 
ohmic A ( 1(0) T(t)) ) 
where 


p = average resistivity of the device calculated from the 


steady-state considerations, 


H 


ee approximate length of the n-region, and 
Ce; = value of the current once the ohmic drop 


is established. 


The value of AV calculated for our case is about -0.1 volts, which is 
well in agreement with our computer result as displayed in Fig. (9.3). 
As the ohmic drop occurs very rapidly, the displacement current in the 
bulk regions is significant only over very short time intervals and is 
replaced by the majority-carrier drift currents. The majority-carrier 
current is the consequence of the presence of a small electric field in 
the bulk regions after the ohmic drop is established. 

As the transient evolves in time, the particle currents become 
significant in the transition region and start extracting the injected 
holes and electrons from the very vicinity of the metallurgical junction. 
The departure of the minority-carriers from the vicinity of the 
metallurgical junction results in a build-up of the electric field in the 
Bie ion region as shown in Fig. (9.4a). Figure (9.4b) displays the 


build-up of electric field throughout the diode and shows the presence 


Beis snall.elactric field in both the p~ and n-regions. Since the 


OW) — 


5 - 


¢ = yc 


si? .(€ po OF 30 editing) Sesgee Soeeaee ae TOR .9OR 
. m4 Z 


= & > } , acs 


lostinn ot ven nabedies ete 2 qemb aie ate tq stan 


Ya, i vitelaaivet sgeisve ~ 
4 ide: ylit pth Ba "bbet9 

i DB Co otenhewiqgs = % 
fT /4ehp eft to salav * (;701 


ahittgges el 


teeo ino 31g? Bebe Iwasa VS to 9m 


Ay 


Si 2 titywes smo @32v Dason 
“nl Ive? CER esimce cB 3 

yino Jsopoliingia 22 ant 
(list ssleves=eoidetee of3 is bo aE! 


p 


>scf | PED % COMevpsangs ada ak 


air 
tedta aaolgat atud 


— 


a 


enoied aida sito efsideeq s j a-vievs Saotenk33 oi) BA. 5 oem 

ae tind 

beers prt ott: tJ s nalget aptviwesa? add tf Japor: 

datsudet Tantate! vido te viiaboty -hae gee anos3oale ‘lana 
a7 i" 

ef) Fo vith i) ated sges+—o eee aaa Te pas 

aff at Biel) atrisals of) wi~biiuwd: & at kannaaie wotjonit Leobguan 

sie eyetghih Cdp.0) og teed, eer at mnieule ba aorges notitnaias') 


a 
ar 


AOMSEOAG and yore by Sloth, anaed 1th Segttey 858 fel whe abispale to qobhted | 
SI ssgnte 6 Bti9! j Dim, ~ toes dagu al ‘Btert oteioate Late, ek: 
{ : oe : a ; 7 as! 7 an a Pe. hl : 


Ta Li 
t 
ib 


Lo . 


os 


“DIODE No. 1 
R= 109, 


t!=0.50 x1077 
t2=0.19 x1073 
1320.30 x1073 
t4=0.41 x1073 
(5= 0.52 x1073 
15=0.69x1073 


=210 


‘© -30 
x 
XQ 
amd 
LL 
ie «tO 
VY 
Oo 
— 
y 
ae) 0) 
LL) 
-60 
-70 
250 pi ie Wisco | 
0 2.4 4.8 7.2 96 12 


DISTANCE (microns) 


Fig. 9.4a. Evolution of the electric field in the vicinity of the junction 


during the initial period of the turn-off transient, for diode No. 1. 
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Fig. 9.4b Evolution of the electric field throughout diode No, 1 during 


the initial period of the reverse transient. 
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density of the majority-carriers is much greater than the density of the 
minority-carriers, this small electric field gives rise to a very large 
majority-carrier drift current which almost accounts for the total 
current. The electron and hole current distributions throughout the 
diode for one particular time instant are plotted in Pig i X95) tithe 
fact that the majority-carrier current cannot be ignored in our case, 
renders the validity of the analytical models adopted for highly 
asymmetrical p-n eons ao doubtful when applied to heavily doped, 
doubly-injecting diodes like the one under consideration. 

Figure (9.6) displays the overall transient response for diode 
No. 1. The computations were discontinued after 0.2 microseconds, since 
we are only interested in times corresponding to the periodic times at 
microwave frequencies. After an initial drop, the terminal voltage and 
terminal current assume constant values. The spatial distributions of 
the hole density during the transient period are plotted in Fig. Coe 7D. 
Figures (9.8a) and (9.8b) display the spatial distributions of the 
electron density during the transient. The distribution of the carriers 
show the effect of electric field in the transition region which tends 
to accumulate carriers at the edge of the transition region in the quasi- 
neutral bulk regions, and deplete the region close to the external 
It may be stated that this effect tends to lengthen the time 


contacts. 


period in which the terminal current remains constant (constant current 


phase). 


9.3 Results for Diode No. 2 


The parameters chosen for the switching circuit of Fig. (9.1) 


for diode No. 2 are as follows: 
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=3 : 
particular time instant (= 0.55 X 10 ~) during the reverse transient, 


for 


diode No. 1. 
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Fig. 9.7. Hole density distribution at various instants of time during 


for diode No. l. 
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Fig. 9.8a. Electron density distribution at various instants of time 


during the reverse transient, for diode No. l. 
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during the reverse transient, for diode No. l. 
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Vp = -5 volts 
Ve = 0.8 volts 
R= 10 Q 


The initial part of the transient is marked by the build-up of the ohmic 
, a ay ; , 

drop in the p — and n -regions and the center high-resistivity region. 

The build-up of the ohmic drop in all the three regions is governed by 

their respective dielectric relaxation times. The estimated average 

dielectric relaxation times, computed on the basis of charge distribution 


before the reverse biasing potential is applied are 6.25 X Oy fe sec. 


a Sec. and 3. Xx eee sec. for the ee n- and p =repiong, 


eeeRX 10 
respectively. The initial part of the transient behaviour is depicted 
in Figs. (9.9) and (9.10). Figure (9.9) displays the terminal voltage 
and the terminal current for the initial part of the transient. Figure 
(9.10) shows the build-up of the electric field throughout the interior 
of the device. The hole and electron density distributions which are 
Potmsiconiticantly affected in a timevof the order of 1.1 xX i046 Sec.s 
are not plotted. The high reverse current caused by the applied reverse 
bias potential and limited by the external series resistance, is 
initially contributed by the displacement current. The build-up of the 
electric field gives rise to large hole and electron drift currents. 
These hole and electron currents become the major contributor to the 
total current within a time of the order of Or 25a Gs sec. «Such, hich 
particle drift currents start depleting the stored charge in the > et 1 
space-charge region. The depletion of the charge from the pr - n space- 


charge region reinforces the build-up of the electric field, which in 


turn causes the build-up of the terminal voltage and consequent decrease 
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Fig. 9.10. Electric field distribution at various instants of time 


(values indicated in Fig. 9.9) during the initial period of the 


reverse transient, for diode No. 2. 
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in the reverse terminal current. The terminal current decreases to 
about 75% of its initial value in approximately 5.5 Xx fry sec. The 
decrease in the terminal current causes the built-up electric field in 
the —_ n- and nl eeeione to decrease slightly. However, at the pe in 
junction, due to the depletion of the charge, the electric field 
continues to build up slowly. Figure (9.11) shows the overall response 
of the device. The reverse current and, therefore, the terminal voltage 
exhibit a relatively slower rate of change after the initial large 
build-up. The terminal current decreases to 25% of its initial value 
within 8 X ee sec, (Figures (9.12), (9.13) and (9.14) display the 
electric field, hole density and electron density distributions for 
various instants of time as indicated. Figure. (9.12) shows that the 
negative electric field at the af -~ n junction is the main contributor to 
the terminal voltage build-up. At the n - ne junction, the electric 
field first rises and then starts decreasing slowly. Figures (9.13) and 
(9.14) show increasing charge depletion in the vicinity of the aa -n 
Gimetion, but exhibit a slight accumulation in the vicinity of the n — n 
junction. This behaviour of the hole and electron density is rather 
snected. It may be explained in the following manner. The present 


model considers the hole and electron mobility to be a strong function 


of the distance. The carrier mobility shows a large variation in the 


vicinity of the regions of different doping (see Fics. (.1g4a), C7), 


@7-13c) and) (7.13d)). Moreover, the electric field in the vicinity of 


the regions of different doping is greater than in the rest of the 


device. The electric field in the n-region is higher than in the n - 


and a eee due to the ohmic drop caused by the large reverse current 


through the device. As the transient proceeds, the reverse current 
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decreases and so does the negative electric field in the n-region. The 
Purrent=-in the vicinity of then - a junction consists of both the 
hole and electron currents. The displacement current is relatively a 
small component and, therefore, can be ignored for the argument here. 
In the eeeee ion: current is mainly contributed by the majority carriers 
and the minority carrier current is negligible. [In the n-region, 
current is made up of both the hole and electron current components. 

The electron current component being almost twice the hole current 
component. Figure (9.15) displays the distribution of both the hole and 
electron current components in and around the n-region for a particular 
time instant. Now let us examine what happens in the vicinity of the 
n- ~~ junction. The holes drifting towards the n-region under the 
influence of the small negative electric field in the Waa find a 
region of greater negative electric field in the vicinity of the n - n 
junction, though the mobility has not changed. This region of greater 
electric field causes the hole current to increase due to a larger drift 
component. As the holes enter the edge of the n-region, the mobility 
increases and as a result the hole current magnitude further increases. 
When the hole reaches that part of the n-region where the electric field 
once again has a low value, the hole current decreases. This results in 
the accumulation of holes in the vicinity of the n - ne junction. As a 
consequence, a gradient of the hole concentrations is formed at the edre 
of the n — a junction so as to support the terminal reverse current. 


Now let us consider the electrons drifting towards the n - a 


ae : 
junction. As they reach the vicinity of the n-n junction where large 


negative electric fields are present, the electron current increases due 


to a larger drift component caused by the electric field. When the 
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TIME = 0.29 1074 
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transient, for diode No. 2. 


Spatial distribution of hole and electron currents for a 
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+- 
electrons cross to the n -region part of the transition region, the 
mobility decreases, resulting in a decrease of the electron drift current, 
As the electrons travel a little further, they encounter a smaller 
eveetric field which results in making the electron drift current still 
smaller, The abrupt decrease in the mobility at the n - oe interface, 
thus , indirectly aids the electrons and the holes to accumulate in 
, ay , 
the neighbourhood of the n ~n junction. The accumulation of both the 
i e ° e ° cL 4 1 
holes and electrons in the vicinity of the n-n interface results in a 
slight decrease in the electric field in the space-charge region at the 
Tex, ; ‘ ’ Pe 
n-n interface. It may be mentioned that if the carrier mobilities 
are assumed constant throughout the device, the small charge accumulation 
+e x ; ac 
at the n-n interface would not be so obvious. This is one of. the 
(Win S ’ : 
reasons that Benda et al. did not observe any charge accumulation in 
their analysis though they observed that the space-charge region at the 
>’ - n junction was mainly responsible for the voltage build-up across 
: Eee ; 
the device, and the space-charge region at the n-n junction 


Contributed. little. 


9 .4 Conclusions 


In this chapter, the "turn-off" transient behaviour of diode 
No. 1 and diode No. 2 has been presented. The value of the current- 


limiting external series resistance has been chosen so as to enable the 


reverse recovery current to be much greater than the initial forward 


current In the case of diode No. 1, the computations were terminated 


a0, 2° X fo7P sec.. since our interest is confined to times corresponding 
2 dee 


to the periodic times at microwave frequencies. 
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The reverse terminal current in diode No. 1 assumes a constant 
value after the initial decay due to the build-up of the ohmic drop. 
The presence of the electric field in the neighbourhood of the p-n 
junction helps maintain the 'constant current phase’, Moreover, the 
small electric field (positive) also 'pushes' the injected minority 
carriers from the regions close to the external contacts towards the 
Space-charge region. This small positive electric field also gives rise to 
large majority-carrier drift currents in the bulk regions and these drift 
currents constitute the major part of the total current. Therefore, the 
majority carrier currents can no longer be ignored, as was done in the. 
singly-injected diodes case. This effect that the minority-carrier 
current constitutes only a small portion of the total current, slows 
down the extraction of stored minority carriers compared to the case of 
the singly-injected diode. 

In the case of diode No. 2, the space-charge region in the 
wrcinity of the >. - n junction plays the dominant role and behaves as a 
normal p-n junction under 'turn-off' transient conditions. The n-n 
junction space-charge region has little effect. In the initial stages 
of the transient, it helps the terminal voltage to build-up, but at 
later times it shows a slight decrease in the electric field resulting 
from a small accumulation of the mobile-carriers. This unexpected 
behaviour may be attributed to the fact that the carrier mobility changes 


abruptly at the n - nt junction and there is an appreciable electric 


= Ra : 
field present in the neighbourhood of the n - n junction. 
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CHAPTER X 
LARGE-SIGNAL, SINUSOIDAL EXCITATION RESPONSE 


Ho thes chapter, the computed results: for the interasteand 
external behaviour of diode No. 1 and diode No. 2 in a broadband, 


harmonic generating circuit are presented and briefly discussed. 


10.1 Introduction 

The broadband circuit of Fig. (10.1) is one of the simplest 
circuit configurations which can be employed for harmonic eeneracionn | 
The internal and external behaviour of the diodes in such a broadband 
circuit is analysed here. This circuit configuration is chosen to 
minimize the complexities of the external circuit and its interaction 
with the device. The analysis is based upon the numerical solution of 
the physical and mathematical model as discussed in Chapter III and 
Chapter IV. The details of the numerical solution are discussed in 
Chapter VI. The signal frequency is chosen to be 5 GHz. to reduce the 


; : tcro Ere ie sen at 
computation time. However, lower microwave frequencies can be cho a 


the expense of more computer time. The d.c. bias is assumed to be 


equal to 5 volts and the peak to peak value of the applied signal is 


assumed to be 12 volts. The load and source resistances (Ro and RO) are 


assumed equal to 3.7 ohms. This value of resistance is within a factor 


of 5 of the computed series resistance of diode No. 1 at zero volts bias, 


(definition given in Chapter V). 


Various quantities, such as electron and hole densities, 


electric field and diode terminal current and voltage are displayed in 
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normalised units, unless otherwise specified (see Table 4.1). Fourier 
analysis of the current and voltage waveforms yield the fundamental 
frequency contents Ena thus enables us to compute the large-signal 


impedance of the diodes at the generator frequency. 


10,2 Results for Diode No. 1 

The numerical solution for the diode behaviour in a broadband 
harmonic generating circuit of Fig. (10.1) is obtained here by applying 
both the reverse d.c. bias, and the a.c. signal simultaneously to the 
diode which is initially in thermal equilibrium. The results are 
Seeained for the first five a.c. cycles. The terminal current, the 
terminal voltage and the applied signal as a function of time are 
meeplayed in Fig. (10.3). It fs evident from Fig. (10.3) that The 
terminal response becomes almost repeatable in the Line domain after the 
initial transient period (=1.1 X Ores sec.). Detailed snap shots of 
the electric field, and the hole and electron density inside the diode 
during the fourth cycle are displayed in Figs. (10.4a), (10.4b), (10.4c) 


fd) (10.4d).\, Figure (10.5) gives the hole and electron densities for 


various Pane anstants during the fourth cyele. Points to notesare the 


following. 


Gb) The electric field in the very vicinity of the metallurgical 


junction (1 to 2 microns on each side of the junction) controls the time- 


dependent behaviour of the device in this case. The electric field away 


from the junction (about 3 microns and more on each side of the junction) 


is small in magnitude; it thus contributes little to the time-dependent 


behaviour. The high electric field at the junction spreads into the 


bulk when the net applied voltage is negative and contracts back towards 
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Fig. 10.1. Circuit model used for the study of the large-signal, 


sinusoidal excitation response of the diodes. 
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Fig. 10.2. Applied waveform in the-cireult-of Figs" G0.1). 
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the junction when the net applied voltage is positive. 


(2) The bulk regions exhibit large majority-carrier drift 
currents which account for almost the total current flowing. These 
current components result from the small electric field which is built 
in the bulk regions due to the flow of the large current. The 
displacement current is very small throughout the diode except in the 
‘very vicinity of the junction, where it contributes significantly to 


Ene total current. 


(3) No significant recovery mechanism is exhibited for 
retrieving the injected minority~carriers which have moved relatively 
far away from the junction (2 microns or more). Large carrier 
concentration gradients here give rise to a forward particle current, 
The negative electric field here is small and does not give rise to a 
large reverse current, so that the minority-carrier currents have a net 
forward component. Thus, the carriers which are injected beyond the 
feeion of the tiigh electric field tend to travel further towards the 
contacts rather than being retrieved during the reverse part of the 
bea However, the carriers present in approximately 1.5 microns width 


on each side of the junction are effectively retrieved during the 


reverse part of the cycle. 


(4) The large-signal impedance of the diode is an effective 


value at the generator frequency obtained by the Fourier analysis of the 


aes 
terminal current and voltage waveforms. The value of the impedance 


obtained in the present case is 0.971 + j0.51 ohms, This value of 


impedance, however, is not unique and depends upon the operating 
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conditions and the surrounding passive circuitry, 


H).o4Results cor Diode,No. 2 

The terminal current and voltage of the diode as a function of 
time are displayed in Fig. (10.6). Detailed results ate computed for 
only one cycle after the initial transient vanishes because of the 
lengthy computer runs required to compute the results for a few more 
Bee cies, The results obtained for diode No. 1 show that the terminal 
response is almost repeatable mid the initial transient has venta 
Detailed snap shots of the electric field, electron density and hole 
density as a function of distance during this a.c. cycle are displayed 
ee 10) 7a), (10.75), (10.42), (00s 7d), (10.76) Mma COE) eeeure 
(10.8) gives the terminal voltage and current as a function of time 


during the initial transient period. Points to note are the following. 


(1) The electric field exhibits a large negative value in the 
n-region during the reverse part of the net applied voltage. However, 
during the foruatd part of the net applied voltage, an appreciable 
positive electric field is built-up in the n-region. Both the positive 
a negative electric fields built-up in the n-region during the a.c. 


cycle may be aerrabuted 60 the ohmic voltage drop caused by the large 


current flow. 


(2) The injected minority-carriers are kept well within the 


center n—region. The large negative electric field built-up in the 


n-region during the reverse part of the net applied voltage enables 


almost all of the injected minority-carriers to be retrieved. 
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(3) The bulk regions exhibit large majority-carrier drift 
currents which account for almost the total current flowing. These 
current components result from the small electric field which is built 
in the bulk regions due to the flow of a large current. The displace- 

2 “ ; ae a : 
ment current is relatively small in the p - and n -regions. In the 
n-region, the displacement current accounts for most of the total 


current and the particle currents are relatively small. 


(4) The large-signal impedance of the diode at the generator 
frequency as defined in Section (10.2) is computed to be 46.102 + 4374.96 
ohms. It may be stated again that this value of the impedance is not 
unique and depends upon the operating conditions and the surrounding 


Passive circuitry. 


10.4 Conclusions 


The internal and external behaviour of diode No. 1 and diode 
No. 2 in a broadband harmonic generating circuit has been analysed. In 
the case of diode No. 1, no significant recovery mechanism is exhibited 


for the injected minority-carriers which have moved relatively far away 


from the junction (2 microns or more). In the case of diode No, 2 the 


injected minority-carriers are kept well within the center n-region. 
The large negative electric field present in the n-region during the 


reverse part of the net applied voltage, enables almost all of the 


injected minority carriers to be retrieved. The large-signal impedance 


of both diode No. 1 and diode No. 2 at the generator frequency shows 


an inductive character. This result indicates that the classical model 


of a varactor diode consisting of a series combination of a variable 
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CHAPTER XI 


MICROWAVE MEASUREMENTS OF THE STATIC CHARACTERISTICS 


OF A MICROWAVE VARACTOR 


In this chapter, various techniques used to measure the static 
properties of a varactor, namely the series resistance and junction 
capacitance and hence static Q, are discussed. Measurements based upon 
these techniques are presented and compared. A commercially available 


step-recovery diode (Hewlett Packard No. 5082-0251) has been used. 


Pie. Introduction 

For most varactor applications at microwave frequencies, it is 
not adequate to rely on low frequency measurements of the varactors 
parameters; even parasitics may be frequency dependent. Besides, 
as of the series resistance at low frequencies is difficult and 
inaccurate. Microwave measurements are therefore necessary for varactor 


characterization. There are two basic methods of measurement presently 


applicable to waveguide-mounted varactors: 


Soo, 
(1) the transmission loss method due to DeLoach ~ with 


: 20 
variations by Roberts and Wilson °; and 


fee io) Ur ae ajeteed 
(2) relative impedance method due to Houlding with variations 


61 A 62 
by Harrison and Mavadatt ©. 


The transmission method 8,29 requires mounting the diode in reduced- 


1 i ing nsmitted power near the series- 
height waveguide and measuring the tra p 


resonant frequency. The method consists of either holding the bias 


fixed and varying the ereaueneyes or holding the frequency fixed and 


i TRAAY HvAORDIM A a 


(4) booe @oopbiedess APE Tea qeris ‘ok 
oy cojvae dats glewia .setosmet me ote 
copsian  jyecrudetlt eee) sheath seen eee 
S%otved bavegwes hem baseneas sis £ 


L4)~SH0P '. 0m devtolag® S3steal) shoth vs 


Ba: 
ia 36 stot eethate wostav Jew. a 


/ OO Ce walla leat val ne yier ba: a 
coum’ ed Yaa eolsidmeg nove i 
sh wok Js sunssetesr astuoa sift to 
an 0 1zradd dee edemstuansm ovewor she 
bwiged ategd pam, Sh piaet »rotiax t 
igiolsewsy ian aude sa bi: 

ionmiot of 36h hediga sol PRN ila On. 
ban {st soe Lil brs astdok t6 eno} is 


helene (gin pen oat eu ed aan somabagat ovtsalss hie oH . 
: al 3s sahavall naw 1 
+hssuhsdi wh sbOgh of): gen ee Ps ioiisent nate ¢ Hi 

22 aiid! inn, aval hieieeaaalthy a) aninennae Bre sagen gt 
ent’ off gutbbDort! sortlhe he agelass Belden oat ‘panei % 
Hil’ howe? Set) ise ont HEL Lotte: 9 cmmnpett ods. | : hmm . 3 2 
9 ; 5 ‘| > i ha Ad 7 

: - s in a 


223 
, BZ! : 
mecying the bias... Due to the necessity of mounting the diode in 
waveguide, this method is limited in the range of junction capacitance 
and package inductance which may be resonated, 
: P 60,61 

The relative impedance method ~’ ae requires matching the 
diode to a slotted line with a lossless tunable transformer. By changing 
the bias voltage and measuring the VSWR, the diode quality factor may be 
calculated by assuming a series resistance which does not vary with bias 
voltage. However, while measuring high quality varactors, the line 
loss and holder loss cause a large measurement error and correction for 

ay 5 63 
these additional losses is necessary. Sard ~ gave a new procedure to 
correct for these losses. However, this procedure is rather tedious and 
complex. 

Another technique is presented here which considers the diode 
as a dielectric post mounted in the center of the waveguide. Various 

2 ‘ hoes : Paar q 

elements of the equivalent circuit are determined by impedance 
measurements. These measurements are corrected to transform the 
impedance measured in a regular height slotted section to that in the 


reduced-height waveguide, in which the diode is mounted, using the 


: , 65 
procedure given by Oliner . 


11.2 Dielectric Post Approach 


a Reduced-Height Waveguide 


11.2.1 Equivalent Circuit of the Diode in a Reduced-Height Waveguide 


An account of the theory pertaining to the situation where the 


coaxial line presents an arbitrary impedance at the junction with the 


waveguide is given by fee ae as. shown..inBig..Glb.l).. Yor simplicity, 


x i Z, 
the coaxial line in Fig. (tl eis replaced by a lumped impedance 1 


connected to the end of the perfectly conducting post, as shown in 
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rectangular 
waveguide 


Ss 


coaxial line 


Fig. 11.1. Sectional view of junction between a rectangular waveguide 


(TE1Q mode) and coaxial line, 


Fig. 11.2. Equivalent case of configuration shown in Fig. CHeTTo)ge 


waveguide 


aye 


Fig. 11.3. Equivalent circuit for the configuration shown in Fig. (11.2). 
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Pee. Cli. 2)%.When. a TE | g7mode electromagnetic wave is incident on a 
perfectly conducting post mounted in a waveguide, the incident electric 
field is cancelled at the surface of the post. Surface currents are 

set up on the post, and their radiation into the waveguide is responsible 
for the reflected wave that is set up. For a post of small diameter, 
tyat is, << 1, where d is the diameter of the post, the resulting 
effect is similar to that produced by a fictitious current filament 
located at the center of the res A solution for the integral 
equation of current yields a value for the reflected wave which is 


6 
dependent on and related to the impedance of the post as follows 


1X l 1 ~1 
Z=2Z cosec* —— cme 1 eg = eis pee Ct Td) 
e a Z j8 p 
1 d 
where 
jB, = capacitance susceptance term, and 
jX_ = inductive reactance term. 
P 


In the course of actual numerical calculations, a value for the applicable 


waveguide characteristic impedance zZ resulted in “: 
r 

= 240 1 {8 eeGulee) 
s a\A 


guide wave-length, and 


N 
| 


where 


> 
ll 


& 


} = free space wave-length. 


int ted b aga) as the 
° - s e 
The capacitance contributing to jB, was interpre y 


capacitance between the post and the walls of the waveguide as well as 


being due to the discontinuity at the junction of the coaxial line with 
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the waveguide. The post inductance for a perfectly conducting short~ 


@-rcuited post is given oes 


fae aye 2 rethatt 
i = 3h log (ig aim #2) S tee Gigs) 


where 


S = correction term, and 


Au 
II 


diameter of the post. 


Figure (11.3) gives the circuit representation of Fq. (11.1). 

Similarly, the impedance for the case of a dielectric post 
instead of a conducting post is obtained by replacing the electric field 
due to the filament current by an appropriate expression for the electric 


field. For a small diameter dielectric post, the post reactance si is 


: 66,67 
given by : 
2 
A WI, (md/A) 
fo) §= cosec*—— fee pre eae eee ee ae 
: er ea ifmavet/ry pau 
nb sj {rave"/ ue J, (md/a) 
paper (SETS 
where 


J_ and Jy = Bessel functions, 
fo) 


' = relative permittivity of the post material, 


E = 
= permeability of the post material, 


w = frequency in radians per second, and 


S = correction factor. 
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Usine—the- explanation: of Getsinger’’. the terms on the right side of 
Eq. (11.4) can be represented by a geometric reactance oS of the post 
in the waveguide, the radial-line admittance of the dielectric post and 
the negative value of the radial line admittance we, at the post 
radius, but with the post absent. If the admittance Y, is complex, 

is replaced by a impedance Zo | 


= apes ee ee See 
Z cosec® —— jwL + ui baa 


x. = jue, 

These considerations lead to an equivalent circuit for a dielectric post 
mounted in a waveguide and terminated in a impedance Zy (1.¢@,,4- coaxial 
cable, etc.). Figure (11.4) gives the complete equivalent circuit. At 
a given frequency, the real and imaginary parts of the post dielectric 
constant can be adjusted (conceptually at least) to provide any desired 
value of radial line admittance over the surface of the post. Let the 
dielectric post be replaced by a diode in a package consisting in part 
of an outer dielectric sleeve having the same diameter as the post and 
having the same radial line admittance over its surface. It would then 
be expected that the waveguide equivalent circuit for the dielectric 


post can also be used to represent the impedance of this new 


configuration with reasonable accuracy. If the diode is centrally 


mounted, which is the usual case, the transformer reduces to one with a 
> 


unity turns ratio. The capacitance C, 1s given by : 
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waveguide 
terminals 


LZ 


7X) 
cosec ~—"] 


ne) 
Fig, V4. Complete equivalent circuit for waveguide shunted by a 
small-diameter dielectric post and an impedance Z, presented by a 


@eoaxtal line; 


rectangular 
waveguide 


Le 


Fig. 11.5. Simplified equivalent circuit applicable for a centrally- 


mounted diode in the waveguide. 
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For a usual microwave diode, Ca is of the order of 0.02 pF. The 
presence of this parameter is ignored for our case as the capacitance 
to be measured is at least 20 times greater than C.. Thus, the 


d 


equivalent circuit reduces to a simplified circuit as shown in figure 


ES). 


11.2.2 Measurements 

To determine the packaged diode impedance from the equivalent 
circuit shown in Fig. (11.5), impedance measurements were made on a 
varactor diode mounted in a reduced—-height rectangular waveguide whose 
height is chosen to be equal to the ceramic portion of the diode. Since 
no slotted line was available in reduced-height waveguide, measurements 
were made in a regular height waveguide (RGU/50) slotted line. The 
reduced—height diode mount was connected to the regular height waveguide 
by an impedance transforming section. The transformer used here was a 
quarter-wave stepped impedance transformer, which was theoretically 
designed to exhibit a return loss greater than 30 db all over the band. 
The actual transformer exhibited a return loss greater than 27 db over 
the frequency band of interest (6 - 8 GHz.) and 29 db for (6.35 — 6.9 GHz) 
along with 30 db at 6.4 and 6.7 GHz + 0.1 GHz, when two identical 
transformers were put back to back. 

The reduced height diode mount is shown in Fig. CLIS iy eee 


uses a coaxial capacitor (more than 100 pF) to serve as a d.c. block 


ppd ar. 6. by-pass Steps were taken to ensure that the contacts at 


both ends of the diode package were practically loselcee, = 1heg. as 


washers were used to allow for small variations in the diodes dimensions. 
A o 


Both the diode mount and the impedance transformer were gold-plated to 
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Fig. 11.6. Cross-section of a typical packaged varactor diode. 


DC BLOCKING 
CAPACITOR 


FLANGE 
‘ 


CENTER SECTION CAP 


REDUCED HEIGHT 
. WAVEGUIDE 
pres 

© 


LM eezeere IZ 
eee ZZ ZT Ze 


Fig. 11.7. The reduced-height rectangular waveguide diode mount. 
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minimize the losses, 

First, the measuring system was ePiGiracea to yield corrections 
for the discontinuity effects introduced by the impedance transformer, 
since the impedance measurements were made with a regular height 
slotted line. A procedure outlined by Ol ineroe allows correction of the 
measured data by means of a calibration curve. To obtain this 
calibration curve, positions of the voltage standing wave minima in the 
miocted line are accurately measured for various positions of an 
adjustable short circuit in the reduced—-height section. These positions 
are closely spaced and readings must be taken over a distance of at 
least one-half of a guide wave-length at the frequency of the r.f. 
signal. From these measurements, a curve showing the difference between 
the position of the short circuit and the position of a minimum in the 
slotted line (both with arbitrary reference planes) as a function of the 
position of the slotted section is Borarneds The calibration curve is 
strictly valid only for lossless terminations. However, Oper © gave a 
procedure in which complex impedances can also be dealt with, using the 


same calibration curve and certain approximations. The applied equations 


then pacomeees 
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where 


D = position of voltage minimum, 


AD = correction available from calibration curve, 
es measured VSWR, 
ae = guide wave-length of moveable short, 
Na = guide wave-length of slotted line section, 
g 
D = position of voltage minimum after correction, 
corrected 
cn = corrected value of VSWR, and 
corrected 
AD = correction available from calibration curve for 


1/8 
a position D - do 


A schematic of the experimental arrangement is shown in Figure 
(11.8). The experimental procedure was as follows: The position of the 
voltage minimum was noted for many closely-spaced positions of the 
adjustable short circuit in the reduced-height section. Readings were 
taken over a distance of more than one-half of a guide wave-length, at 
the frequency at which measurements were desired, The adjustable short 


circuit was placed at a distance Pl behind the position in which the 


diode is mounted. The diode and coaxial line were then removed and the 


resulting hole in the broad wall of the waveguide was closed by a gold- 


plated flat brass plate. The position of a reference voltage minimum in 


the slotted line was noted and corresponded to a short circuit at the 
plane of the diode. Then the adjustable short circuit was placed at a 


distance A /4 behind the location at which the diode is mounted, thus 


corresponding to an open circuit at the plane of the diode. With this 


position of the short curcuit, the VSWR and position of the minima was 
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noted when the coaxial line was reinserted with a highly conductive post 
of the same shape and size in place of the diode. The conducting post 
‘was then replaced with the diode, and the position of the voltage 
minimum and the VSWR were noted, the short circuit remaining in the 
Same position. These readings were taken for various d.c. bias voltages 
applied to the diode. Ail these measurements were corrected as given 
in aera’. ) and (1c) > From the equivalent circuit siven im fie. 
(11.5), it is seen that the diode terminal impedance at any bias is 
readily calculated from the above measurements as the difference 
between the measured impedance of the assembly and the Memes of the 
system with a aummy diode substituted for the varactor. A calibration 
curve for a frequency of 6.757 GHz is shown in Fig. (11.9). ee Srna 
on the uP step-recovery diode No. 5082-0251 are listed in Table (11.1). 
The value of the Pee eee seas impedance Zo calculated from Eq. (11.2) 
foe 55.0 Olms Lor this’ case. 


25 Obtaining the Series Resistance from the Measured 


Terminal Impedance 


The most commonly accepted equivalent meni of a varactor 
whose cross-section is as shown in Fig. (10.6) is as shown in Fig. (11.10), 
for the small signal case. For the measurements at the frequency 
considered here, fringing capacitance Ca can be neglected and the 


equivalent circuit simplifies to the circuit shown in Fig. (11.11). The 


value of the package inductance Ly and package capacitance Cy can be 


67 : 
measured as outlined by Getsinger . However, to perform the 


measurements of L, and Cy with sufficient accuracy, we need a set of 
A 


a a 


peereer 


: 68 
* Similar technique was applied to tunnel diodes by Scanlan . 
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TABLE 


Corrected 


Relative 


Impedance 


0.030 + 50.482 


0.067 


0.058 


OG 34.5 


O2.053 


KsOrt 


0-079 


ey tea atone aaranuesenabengienssiimasemsasweushmpssieen 


Relative Diode 


Terminal 
Impedance 


Q5037-="70" 

0,028 =. 40. 
0.024 + 70. 
023 = 40, 
Os0a2-= 40). 
0.049 + 0. 
el Oe 


70 24 Oe 
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015 


368 


453 


True Diode 
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Impedance 
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(HP: NO. 5082-0251) MEASURED BY THE 


DIELECTRIC POST APPROACH. 


.1. TERMINAL IMPEDANCE OF STEP-RECOVERY DIODE 
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Fig. 11.10. Packaged varactor diode equivalent circuit at microwave 


frequencies. 
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Fig. ll.11. Simplified equivalent circuit of packaged varactor diodes 


at lower microwave frequencies. 
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open and shorted diode packages so that statistical data can be collected 
from ee a usable value of Cy and Ly can be chosen, Sines we do not 
have a set of open and shorted diode packages, we needed to rely on the 


specifications supplied by the manufacturer for the package: 


OG 
Ht 


0.20 pF 


Le = 0, 50Tan 


From the equivalent circuit of Figs (laa tye we can write the 


following equation 


So st ae ee ee ne ae 2. ALT. 9) 


where 


ve) 
H} 


real part of the diode terminal impedance. 


If it is assumed that the junction capacitance a at the measurement 
frequency has the same value as determined by the low frequency 
* OOGG3) me 


capacitance measurements » Ro can be determined directly from 


Eq. (11.9), since we already know the value of La and Cy. 


Silk3 Results 
The diode terminal capacitance, measured with a Boonton bridge 
using a small 1 MHz signal to minimise averaging effects, is displayed 


myrte. (11,12) as a function of the bias voltage. From the equivalent 
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eeccuit of Fig. (11.11), this value of capacitance is almost equal to 
the value of (Cc, z C,). If it is assumed that the junction capacitance 
oe at the measurement frequency (6.757 GHz) has the same value as 
determined by the 1 MHz frequency measurémentee >>>, the junction 
capacitance is determined from the measurements reported in Fig. (11.12), 
since we know the value of Cy (0.2 pF) as determined in Section (11.2.3). 
Equation (11.9) then yields the value of the series resistance Ros since 
we know all the other parameters of the equivalent circuit. Table (11.2) 


lists the value of the series resistance R. determined from Eq. (11.9) 


for various bias voltages. 


Bias Voltage Series Resistance 
in volts. d.c. 3 Rg in 2 


3 5 TEP-RECOVERY DIODE 
TABLE 11.2. THE VALUES OF SERIES RESISTANCE R, FOR S$ 


(HP. NO. 5082-0251) AS DETERMINED BY THE DIELECTRIC POST 
APPROACH. 


Figures (11.13) and (11.14) give the static Q measurements at 


6 GHz and 7 GHz, respectively, on the same diode using Harrison's 
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Step-recovery diode (HP. 5082-0251) impedance for various bias voltages derived by 


Harrison's method at 6 GHz. 


Fig. 11.13. 
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relative impedance method. A shorted diode package provided by Hewlett 
Packard Company is used to obtain the required short circuit reference, 
This method yields Static Q of 7.7 at 6 GHz and 6.5 at 7 GHz at zero 
volts bias for the diode under consideration. The static Q at zero 
volts bias for the dao aet: calculated by the dielectric post approach 
meeoented in Section (11.2) is equal to 8.756 at the measurement 
frequency of 6.757 GHz. Figure (11.15) compares the value of the 

Pe itance as calculated for diode No. 2 for steady-state cases and the 
values of capacitance obtained for the step-recovery diode (HP. No. 
5082-0251) from the Boonton bridge operating at 1 MHz. The values show 
an agreement for low forward-biases. The relatively large difference 
at high forward-bias may be attributed to the diffusion capacitance 


dependence upon the alee of the measuring signal. 


11.4 Conclusions 

An improved procedure has been developed for measuring the 
static characteristics of a varactor diode. First, a theoretical circuit 
model for a varactor mounted in a reduced-height waveguide has been 
given. Then an experimental procedure for determining the parameters 
of this equivalent circuit has been described. This method takes into 


account the line loss and holder loss which are present in a practical 


measuring system. This method also yields directly the value of the 


Beries resistance R.. The value ofethe static Q measured from this 
Ss 


approach is higher than the one measured by Harrison's approach as the 


1 
latter approach neglected line loss and holder loss, and thus assumes 


the measuring system to be lossless. 
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footie No, 2 


mam MEASURED VALUES FOR 
hp 59082-02511 
STEP-RECOVERY DIODE 


INCREMENTAL CAPACITANCE (p,) 


0 Ty: 0.4 0.6 0.8 1.0 
DIODE TERMINAL VOLTAGE (volts) 


Fig. 11.15. Calculated incremental capacitance of diode No, 2 and low 


frequency terminal capacitance of step-recovery diode 


(HP. No. 5082-0251) as functions of forward-bias voltage, 
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CHAPTER XIL 
CONCLUSION 


In this work, a computer-aided study of the electrical behaviour 
of microwave varactor diodes under static and dynamic large-signal 
conditions has been carried out. Numerical methods of solution have 
been used in preference to presently available analytical techniques 
because of the freedom available in the former approach in the choice of 
doping profile, generation-recombination law, mobility dependencies and 
injection level. Moreover, boundary conditions in this method are 
applied solely at the external contacts. 

For purposes of this study, microwave varactors have been 
classified into two broad categories, namely varactor diodes and step- 
recovery diodes. The first category was taken to represent those diodes 
in which non-linearity of the junction capacitance in the reverse-bias 
condition is optimised. Diodes optimised for utilising the abrupt 
change of capacitance for voltage near the contact potential were taken 
to belong to the second category. Typical physical structures and 
doping profiles have been chosen to represent these two types of 
The steady-state solution has been obtained by an algorithm 


varactor. 


origionally used by Guanelae to solve the steady-state case for 


transistors. Time-dependent solutions have been obtained by finite 


aint : : 
difference methods originally used by DeMari for a simpler physical 


“42 x -?t ° 
model of zero recombination and constant mobility. DeMari's algorithm 


has been extended to cover our physical model and was successfully 
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applied to the two different chosen diode structures; they exhibited 
a very different stability behaviour in their numerical solution, the 
step-recovery diode solution requiring more stringent procedures. 

Four’ different forms of external excitation have been 
considered: 

(1) steady-state voltage excitation; 

(2) large-signal "turn-on" excitation; 

(3) large-signal "turn-off" excitation; 

(4) large-signal sinusoidal excitation in a broadband circuit. 

Detailed results for the above mentioned cases have been 
presented and. discussed. in Chapters VII, VIII, IX and X. The results 
reflect the usefulness of numerical techniques for studying both the 
static and dynamic behaviour of these devices. These results also 
expose various inaccurate assumptions employed in much of the analytical 
work and in some cases, their inapplicability to the present device 
under the conditions studied. 

An improved technique to measure the static characteristics of 
a varactor diode at microwave frequencies is reported in Chapter XI. The 


value of static Q obtained by this "dielectric post approach' is higher 


61 : . ; 
than given by Harrison's method which neglects the circuit losses. 


This investigation can be used as a basis for studies involving 


the device optimisation and optimisation of the device together with its 


embedding circuitry. A few examples are listed below: 
(1) the effect of doping concentration , width and gradient 


of the net doping profile at the boundaries of the center high 


resistivity region on the reverse-transient characteristics of step- 


recovery diodes; 
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(2) the dependence of large-signal impedance at the generator 
frequency upon signal amplitude and on frequency in a broadband CLECULE® 

(3) the simulation of the circuit conditions when only 
certain frequencies are allowed to flow; 

(4) the extension of the present one-dimensional model to two 
or three-dimensional cases to study complex geometry effects on the 
device behaviour; 

(5) measurement of large-signal impedance of the device at 


microwave frequencies. 


Several extensions of the computer model developed here come 


to mind-and are:Llisted below: 


(1) The present numerical scheme may be extended to study the 
f : + a5 
large-signal behaviour of the single-drift region (p -n-n ) 
; + + 
avalanche diode, the double-drift region (p - p-n-n ) avalanche 


diode and microwave transistors; 


(2) the present scheme may be extended to study the effect of 
multi-level carrier trapping and recombination on the behaviour of 


microwave semiconductor junction devices; 


(3) the present scheme may also be extended to study the effect 


of temperature on the steady-state behaviour of microwave semiconductor 


junction devices. 
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APPENDIX A 


AUTOMATICALLY ADJUSTABLE SELECTION OF 


DISCRETIZATION MESH SIZE IN THE TIME DOMAIN 


ros desired to achieve an anesaae ie adjustment of time step 
size subject to the requirement of maintaining constant errors, throughout 
the time domain, in our numerical solution of time-dependent cases as 
reported in Chapter VI. In the procedure adopted in Chapter VI, only 
numerical differentiations are present. The total error introduced by 
such operations may be regarded as a sum of truncation and round-off 
errors. For a decrease in the step size and an increase in the order of 
interpolation scheme used, the former type of error decreases whereas the 
- latter increases. For the present case, round-off errors may be usually 
considered negligible since double precision arithmatic is employed. 
This is a consequence of practical restrictions on both the order of the 
numerical scheme and the minimum step size that may be used (limited by 
memory size and computation time). The effort will then be towards 
achieving a constant truncation error throughout the time domain. 

The Taylor expansion about a point t = we of a continuous 


function f(t) with continuous derivatives is considered: 
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If a two-point formula is used to approximate the first derivative at 


ene point it = wee such as 
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and is inserted in the expansion (A.2) evaluated at the point t = pte 
one obtains: 
1 Z foe 
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If higher order contributions of Eq. (A.3) are neglected, the term 
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represents the error introduced in the approximation (A.2). The relative 


eats 


error is then defined as the ratio ee / gone 
k+1 


repeated for a sequence of adjacent points ,t Cesena 2 aye toe 
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the requirement of maintaining a constant magnitude of the relative error 


at each point may be specified. In this case the following relation must 


hold: 
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which represents a relationship between consecutive steps in terms of 
local increments and curvatures of the function f(t). 

In the case of interest, the magnitude of the time step pp At 
at the next instant of time may be approximately predicted with the 


knowledge of quantities only available for times t = ,t if the expression 


k 
(A.4) is modified to the following form 
eee 
k-2 k 
At = ( aa ae ee At; (k=3,4,5, ° ° =) ee OS). 
k+1 weep ot k 


It is apparent that the sizes of the initial three time steps (At, 
hoie7.3, « 2.) must be independently specified, and the size of the 
fourth (and following ones) may be predicted by the relation of Eq. (A.5). 
Instants of time for which the argument of the square root of Eq. (A.5) 
exceeds the permissible range become singular points. These are likely 
to occur during the evolution of a transient solution and are easily 
taken care of by specifying the upper and lower bounds for the ratio of 
the consecutive steps and for the magnitude of the step size. These, as 
well as the initial three steps, must be chosen with regard to the total 
number of instants of time that may be allowed and to the rate of 


convergence of the iterative procedure, very much dependent upon the 
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step SUze, 

Although a time step adjustment of the described type is not 
very accurate, it appears satisfactory from any practical point of view, 
also in consideration of the insignificant computation load required and 
by the observation that the selection of the time step is by no means 


ericical,. 
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